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The  capabilities  of  a  Swept  Line  Electron  Beam  (SLEB)  in  annealing 
lon-lmplanted  semiconductors  are  examined.  This  technique  employs  a  fixed 
geometry,  line-shaped  electron  beam  through  which  implanted  samples  are 
mechanically  scanned.  In  general,  this  technique  can  produce  annealing  re¬ 
sults  comparable  or  superior  to  those  achievable  by  conventional  furnace 
annealing. 

Residual  point  defects  in  self-implanted  amorphous  silicon  treated  by 

SLEB  and  furnace  processes  are  examined  by  Deep  Level  Transient  Spectroscopy 

Despite  high  temperature  treatment,  furnace  annealed  samples  show  large 
16  -3 

(10  cm  )  defect  concentrations  and  dopant  migration  phenomena.  This  is 
especially  true  in  the  as-lmplanted  amorphous-crystalline  transition  region. 
When  proper  annealing  parameters  are  used,  SLEB  annealed  material  shows 
much  reduced  point  defect  concentrations  and  reduced  dopant  motion.  These 
relatively  thick  amorphous  layers  (0.5  pm)  are  regrown  and  annealed  by  SLEB 
without  the  use  of  additional  furnace  treatment. 

Similar  studies  of  BF^  implanted  silicon  are  also  presented.  Differ¬ 
ential  reslstivlty/Hall  effect  and  Secondary  Ion  Mass  Spectrometry  analysis 
are  used  to  show  improved  electrical  activation  and  only  limited  dopant 
motion  during  SLEB  annealing  as  compared  with  furnace  annealing.  The  im¬ 
proved  electrical  activity  is  especially  significant  in  the  original 
amorphous-crystalline  transition  region  where  reduced  residual  defect  den¬ 


sities  are  observed. 


SLEB  annealing  effectiveness  in  both  direct  and  indirect  band  gap 
composition  GaAs^_x?x  is  also  Investigated.  Photo luminescence  emission 
from  nitrogen  implanted,  beam  annealed  material  is  found  to  be  comparable 
or  larger  in  intensity  when  compared  with  optimally  prepared  furnace 
annealed  material.  Photo luminescence  profiling  and  p-n  junction  studies 
show  that  migration  of  implanted  nitrogen  and  related  damage  can  be  limited 
to  the  as- implanted  profile  when  SLEB  annealing  is  used. 
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1.  INTRODUCTION 


1.1  Overview  of  Beam  Processing  Technology 

The  use  of  concentrated  beams  of  photons,  electrons  or  ions  to 
prepare  semiconductor  material  has  recently  become  an  area  of  intense  re¬ 
search.  Although  lasers  have  provided  other  Industrial  applications  for 
many  years,  beam  techniques  have  only  recently  demonstrated  the  potential  to 
revolutionize  semiconductor  device  fabrication.  In  this  section,  I  outline 
recent  developments  and  applications  of  beam  processing,  point  out  major 
trends  and  describe  the  relationship  of  my  research  to  the  mainstream  of 
technology. 

The  techniques  used  for  beam  processing  are  distinguished  by  their 
beam  source  and  method  by  which  the  beam  energy  is  directed  onto  the  semi¬ 
conductor.  To  date,  six  beam  processing  techniques  have  been  studied. 

These  employ  pulsed  lasers  (PL),  raster  scanned  CW  lasers  (CWL),  pulsed  ion 
beams  (PI) ,  area  pulsed  electron  beams  (PEB) ,  raster  scanned  (RSEB)  or  multi 
scanned  electron  beams  (MSEB) ,  and  swept  line  electron  beams  (SLEB) .  A 
general  overview  of  experimental  processing  parameters  is  presented  in 
Figure  1.1.  Here  total  average  incident  energy  density  is  plotted  against 
process  duration.  In  cases  of  multiple  scanning,  energy  density  is  summed 
after  many  passes  by  the  beam.  Data  here  are  by  no  means  all  inclusive,  but 
reflect  the  general  range  of  each  processing  technique.  Several  trends  are 
noteworthy.  First,  there  is  a  tremendous  span  over  which  successful  beam 
processing  applications  have  occured;  eleven  orders  of  magnitude  in  time, 
seven  orders  in  energy  density.  Much  of  the  initial  experimentation  has 
concentrated  on  pulsed  techniques,  particularly  pulsed  lasers  and  the  SPIRE 
area  PEB  machine.  (See  lower  left  of  Figure  1.1.)  Processing  is  generally 
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performed  In  a  single  pulse,  thus  limiting  experimental  parameters  to  the 
characteristic  laser  or  electron  discharge  decay.  In  more  recent  studies, 
scanned  continuous  lasers  and  electron  beams  (MSEL,  SLEB,  CWL)  expand  the 
range  of  study.  Scanning  speed  and  repetition  rate  enter  as  additional 
parameters  and  extend  incident  energy  densities  1000  fold,  and  process 
duration  100,000  times. 

By  far  the  most  widely  studied  application  for  all  types  of  beam 
processing  is  the  annealing  of  ion-implanted  layers  in  semiconductors.  Other 
applications  include  induced  crystallization  of  amorphous  and  polycrystalline 
layers,  pulsed  diffusion  of  dopants  directly  from  a  surface  layer,  contact 
formation,  growth  and  Q  reduction  in  oxides,  surface  smoothing,  circuit 

S3 

customizing,  trimming  of  thin  films,  scribing  device  chips,  drilling  and 
mesa  formation.  Recently,  beam  processing  techniques  have  improved  to  the 
point  that  device  fabrication  has  been  attempted.  To  date,  resistors, 
capacitors,  diodes,  bipolar  transistors,  MOSFETs,  and  solar  cells  have  been 
fabricated  with  varying  degrees  of  success.  A  representative  sample  of  pro¬ 
cessing  applications  is  presented  in  Table  1.1. 

1. 2  Beam-Semiconductor  Interactions 

Interaction  of  photons  and  electrons  in  solids  is  in  itself  an 
area  of  active  research.  However,  it  is  appropriate  that  we  review  in 
general  terms  the  important  aspects  of  photon  and  electron  energy  deposition 
in  semiconductors  [51-53].  There  are  three  areas  of  concern:  energy  deposi¬ 
tion  depth  profile,  parameters  which  affect  this  profile,  and  beam- induced 
effects. 

Photons  (0.1  -  5.0  eV)  transfer  energy  to  a  crystalline  lattice 
through  Induced  vibrations  of  valence  electrons.  Reflection  arises  from 


TABLE  1.1 


EXAMPLES  OF  BEAM  PROCESSING  APPLICATIONS 


APPLICATION  PROCESS  USED  REMARKS  REFERENCES 

-Annealing  of  Ion- 
Implantation  Induced 
Damage/Dopant 
Activation 

*  Si  Crystal  Substrates 
**Implanted  Species 


In,Al,Ga,P,As,B 

Ruby  PL 

1-4 

Ga,Ge,As,Se,In,Sb 

Nd : Tag  PL 

5 

Sn,Te,Bi,P,As 

CO2  CWL 

MeV  Ion-Implants 

6 

As 

Argon  CWL 

EBIC. Damage  Invest. 

7,8 

B 

Krypton  CWL 

8 

P,As,B 

MSEB 

9-11 

*  Si  Polycrystal 

**Dopant 

B 

Dye  PL 

12 

P 

Ruby  PL 

13 

*  GaAs  Crystal  Substrates 

**Implanted  Species 

Te,SiySe,Kr,Zn,Mg 

Ruby  PL 

14-21 

SlyTe,Se,KryBe 

SPIRE  PEB 

16-18, 

-Recrystallization 

* 

of  Silicon 

Ruby  PL 

Deposited  Si  on  Xtal 

24 

ND:Yag  PL 

Deposited  Si  on  Xtal 

&  Poly-xtal  Subst. 

25-27 

Argon  CWL 

Deposited  Si  on  Xtaly 

■ 

Lateral  Epitaxy 

8,28-: 

Spire  PEB 

Deposited  Si  on  Xtal 

24 

SLEB 

Ion-Implanted  Si 

31 

-Direct  Doping  from  a 

Surface  Layer 

*S1  Substrates 

CO,  CWL 

Photo-Induced  Zone 

4 

Migration 

32 

*GaAs  Substrates 

SPIRE  PEB 

From  AsSe  Layers 

33 

I 


5 


1 
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TABLE  1.1  (Continued) 


if 

» 

k  *,» 

APPLICATION 

PROCESS  USED 

REMARKS 

REFEREN 

•  ;•;! 

«  ,• 
i 

-Contact  Preparation 

! 

*In  Silicon 

Ruby  PL 

Si/CoSi2/Si  Contacts 

34 

»•  .  ^ 

»  •*” 

ArgonCWL 

8 

•'  *  .* 

PEB,  MSEB 

A1  Contacts 

35 

* 

*In  GaAs 

Nd:Yag, 

,«  '.T* 

Ruby  PL 

Au:Ge  Contacts 

36 

|  * 

Argon  CWL 

In-Au:Ge  Contacts 

37 

*  »'/. 

4  %’* 

-Beam  Induced  Oxide 

•  * 

Growth 

* 

*Si:Si02,Si3NA 

Argon  CWL 

0,N,C  Implanted  Si 

38 

L  V 

1  E 

*GaAs:Ga203 

Argon  CWL 

39 

" 

-Oxide  Charge  Reduction 

:  »*  i 

; 

« 

In  Si02 

Argon  CWL 

8 

!  1 

-Device  Fabrication 

-r* 

• 

*Si  Diodes 

> 

/  rw 

Ruby  PL 

10 

Nd:Yag  PL 

40 

Argon  CWL 

40,41 

1 

MSEB 

10,42 

1  s 

FEB 

40 

;  w% 

< 

*Si  Bipolar 

• 

Transistors 

Ruby  PL 

B,As,P  Implanted 

46 

! 

*Si  MOSFETs 

Nd:Yag  PL 

On  Poly-Silicon 

27 

i 

Argon  CWL 

Ion- Implanted  on  Xtal 

47 

*  ‘  1- 

; 

•  \  ■  * 

MSEB 

91  II  II 

48 

k  •  ■-* 

» 

» 

*Si  Solar  Cells 

Ruby  PL 

P+  Ion-Implanted 

49 

■  !>: 

i  U. 

Ruby  PL 

Pulsed  Diffused 

50 
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partial  reradiatlon  by  excited  electrons.  The  remaining  energy  is  transfer¬ 
red  to  the  lattice  via  electron-atom  collisions.  For  normal  incidence, 
power  density  absorbed  at  a  depth  z  in  a  one  dimensional  semi-infinite 
sample  illuminated  by  a  beam  of  intensity  IQ  is  given  by: 

*(z)  -  Io(l-R)oe"az  (1.1) 

where  R  is  the  fraction  of  photons  reflected  and  a  the  absorption  coeffi¬ 
cient.  Electrons  (5—50  keV) ,  on  the  other  hand,  may  interact  either 
elastically  with  the  atomic  nuclei,  or  inelastically  with  local  electrons. 

The  elastically  scattered  (backscattered)  electrons  can  reemerge  from  the 
sample  with  as  much  as  50%  of  the  energy  of  the  incident  beam.  Inelastically 
scattered  electrons  are  usually  modeled  as  losing  their  energy  to  the  lattice 
in  a  continuous  scattering  process: 

-dE/ds  «  (n  /E  )  ln(E  /J)  (1.2) 

CO  o 

where  Eq  is  the  Incident  electron  energy,  ng  is  the  total  density  of  elec¬ 
trons  and  J  the  ionization  potential  of  the  material.  This  formulation 
applies  along  the  path  of  each  beam  electron.  Since  each  electron  is  scat¬ 
tered  many  times,  deposited  energy  depth  profiles  must  account  for  deviation 
in  the  electron  trajectory  due  to  each  collision.  The  results  of  these 
calculations  are  truncated  pseudo-gausslan[52,53]  distributions  which  are 
relatively  independent  of  the  beam  energy  when  plotted  as  a  percentage  of 
extrapolated  electron  range  [53].  For  silicon,  this  range  is  given  by: 

R  -  4.57  E  1,75  (1.3) 

g  o 

where  E  is  the  electron  energy  in  keV,  and  the  R  is  given  per  unit  density 
^  8 

2 

of  the  target  material  in  yg/cm  .  The  energy  loss  distributions  become 
more  skewed  towards  the  surface  of  the  sample  as  the  atomic  number  of  the 
target  increases  [52]. 


The  probability  of  photon  absorption  la  related  to  the  density  of 

available  Initial  and  final  states  of  a  valance  electron  transition.  Hence, 

both  photon  absorption  and  reflection  coefficients  depend  strongly  upon  the 

photon  wavelength  and  sample  temperature.  In  crystalline  silicon,  for 

example,  a  has  a  value  of  about  10^  cm  ^  at  wavelengths  shorter  than  0.4  pm 
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then  falls  rapidly  to  a  few  times  10  cm  at  2  pm,  as  the  corresponding 
photon  energy  falls  below  bandgap.  Most  semiconductor  bandgaps  decrease 
with  increases  In  temperature,  causing  the  drop  off  in  a  to  occur  at  longer 
X.  The  degree  to  which  a  sample  may  have  heavy  disorder  or  ion-implantation 
induced  damage  will  also  increase  absorption.  As  indicated  in  Eqns.  1.1 
and  1.2,  electron  power  absorption  is  generally  only  effected  by  material 
density  and  x'  vnic  number,  with  far  less  sensitivity  to  temperature.  In  the 
same  way,  the  backscatter  efficiency  (n)  depends  primarily  upon  the  atomic 
number  of  the  target. 

Figure  1.2  demonstrates  the  relative  absorbed  power  density  pro- 

2 

files  in  crystalline  silicon  for  a  1  J/cm  flux  of  2.3  eV  photons  and  20  keV 

electrons.  These  energies  are  typical  of  beam  processing  experiments.  The 
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photon  profile  is  based  upon  Eqn.  1.1  with  R  -  0.37  and  a  ■  8  x  10  cm 
while  the  electron  profile  uses  n  ■  0.174  and  the  experimentally  verified 
energy  dissipation  distributions  of  Everhart  and  Hoff  [53].  Although  the 
overall  depth  of  penetration  is  nearly  equal,  the  power  dissipated  by 
electrons  exceeds  that  of  photons  for  most  of  the  profile. 

Under  more  intense  laser  illumination,  semiconductor  materials 
display  beam-induced  effects  which  significantly  affect  energy  absorption 
[51].  Free  carriers  generated  by  the  photon  flux  can  become  numerous 
enough  to  Increase  a  by  direct  absorption.  At  higher  energy  densities,  these 


carriers  can  produce  impact  ionized  carriers,  further  Increasing  a.  If  the 
photons  are  delivered  in  a  short  pulse,  electron-hole  densities  can  become 
high  enough  to  create  a  reflecting  plasma.  These  hot  carriers  can  also 
diffuse  from  the  surface,  significantly  affecting  the  absorbed  power  pro¬ 
file.  If  heated  beyond  the  melting  point,  both  a  and  R  increase,  due  to  the 
metallic  nature  of  the  molten  state.  Since  electron  power  absorption  does 
not  depend  upon  the  density  of  electron  states  in  the  target,  analogous  high 
power  effects  are  not  significant  for  electron  beams.  Photons  used  in  some 
beam  processes  have  wavelengths  equal  to  major  fractional  parts  of  the  thick¬ 
ness  of  dielectric  coatings.  The  result  can  be  constructive  interference 
and  severe  localized  heating.  Corresponding  electron-dielectric  effects 
are  not  possible  in  the  scale  of  integrated  circuit  structures. 

1.3  Trends  In  Beam  Processing  Technology 


Semiconductor  beam  processing  technology  is  still  very  much  in  a 
developmental  state.  However,  enough  is  known  about  the  characteristics  of 
each  technique  to  predict  its  ultimate  industrial  applicability.  The  re¬ 
quirements  of  each  application  will  be  the  deciding  factor. 

Activation  of  ion- implanted  dopants  in  single  crystal  semiconduc¬ 
tors,  although  the  most  studied  application,  places  the  most  severe  restric¬ 
tions  on  beam  technology.  Activation  of  the  Implanted  species  must  be 
efficient.  Residual  defect  concentrations  must  be  low.  Conventional  sample 
heating  either  during  or  after  beam  treatment  is  not  desirable.  Although 
most  beam  processing  techniques  have  achieved  activation  of  dopants,  a 
number  of  recent  studies  Indicate  tliat  pulsed  techniques,  in  which  the  semi¬ 
conductor  actually  melts,  produce  residual  defects  which  cannot  be  eliminated 
without  post-processing  furnace  annealing.  One  might  therefore  Imagine  that 


any  fast  melt-refreeze  process  would  leave  residual  disorder  at  the  point 
of  furthest  advance  of  the  melt  front.  This  seems  to  be  borne  out  by  recent 
research. 

Bentlni  et  al.  [54]  have  demonstrated  that  MSEB  solid  phase  re¬ 
growth  In  phosphorus  Implanted  silicon  produces  better  minority  carrier  dif¬ 
fusion  lengths  that  both  PEB  and  melt  phase  ruby  PL  processing.  Anderson 
et  al.  [17]  use  differential  resistivity  techniques  to  demonstrate  consis¬ 
tently  larger  mobility  in  Implanted  GaAs  layers  when  CW  lasers  are  employed. 
Ruby  and  Nd:YAG  PL  and  the  SPIRE  PEB  machine  were  compared  with  argon  and 
Nd:YAG  CWL.  They  suggest  that  rapid  solidification  may  be  accompanied  by 
large  concentrations  of  antlstructural  defects.  They  also  suggest  that 
highly  non-stoichlometrlc  regions  can  be  created  by  melt  phase  Induced 
arsenic  evolution.  Direct  evaluation  of  defects  after  both  melt  and  solid 
phase  processing  has  been  performed  by  Benton  et  al.  [55].  DLTS  studies  on 
B  and  As  Implanted  silicon  show  as  many  as  100  times  (10x  vs.  10  cm  '*) 
more  residual  defects  in  Nd:YAG  PL  as  In CO2  CWL  processed  diodes.  In  a 
separate  study,  Benton  et  al.  [56]  use  DLTS  and  TSCAP  to  identify  defects 
In  single  crystal  material  generated  directly  by  Nd:YAG  PL  Induced  melting. 
Boroffka  et  al.  [10]  have  studied  leakage  currents  In  arsenic  implanted 
silicon  diodes  prepared  by  ruby  PL  and  scanned  e-beam.  They  observe  con¬ 
sistently  smaller  leakage  by  the  solid  phase  e-beam  process.  Likewise, 
Sandow  [40]  reports  superior  I-V  characteristics  in  arsenic  implanted,  ar¬ 
gon  CWL  annealed  diodes  compared  with  both  pulsed  Nd:YAG  and  PEB  processed 
devices.  He  suggests  that  the  differences  are  due  to  large  concentrations 
of  residual  defects  remaining  after  the  pulsed  processing. 
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As  one  might  expect,  melt  phase  processing  is  accompanied  by 
rapid  redistribution  of  dopants  [57,58].  This  characteristic  is  not  desir¬ 
able  in  cases  where  special  profiles  have  been  established  by  multiple  im¬ 
plantation,  or  where  great  control  of  the  junction  depth  is  necessary.  This 
taken  with  the  evidence  about  junction  defects  makes  it  unlikely  that*  pulsed 
beam  technology  will  replace  standard  furnace  diffusion  or  annealing  in 
applications  where  high  quality  or  tailored  profiles  are  required.  Continuous 
wave  lasers  and  multi-scanned  or  swept-line  electron  beams  are  the  most 
likely  candidates  for  this  application. 

In  devices  reouiring  dielectric  coatings,  unintentional  spot  heating 
due  to  laser  interference  effects  can  have  a  devastating  effect  on  surface 
morphology  [59].  These  effects  are  not  expected  for  electron  processes.  In 
addition,  raster  scanned  CWL  techniques  have  been  shown  to  produce  nonuni¬ 
formities  in  annealing  quality  on  the  scale  of  the  beam  spot.  Mizuta  et  al. 
[60]  have  studied  As+  implanted  CWL  annealed  silicon  by  electron  beam  induced 
current  (EBIC)  analysis.  They  observe  reduced  minority  carrier  charge  col¬ 
lection  efficiency  in  banded  patterns  parallel  to  the  direction  of  travel  of 
the  beam  spot.  These  bands  are  attributed  to  laser  induced  defects  formed 
in  uneven  patterns  beneath  the  implanted  layer.  Since  spot  sizes  were  of  the 
order  of  microns,  such  banding  can  have  serious  effects  on  the  characteristics 
of  planar  devices.  Recently  Sheng  et  al.  [7]  have  extended  this  EBIC  study 
to  include  raster  scanned  electron  beam  annealing.  They  demonstrate  that  the 
RSEB  technique  produces  even  less  residual  damage  than  argon  CWL  treatment, 
without  the  banded  damage  regions. 

The  facts  outlined  above  seem  to  Indicate  the  general  suitability 
of  electron  beams  for  processing  semiconductor  material  of  integrated  circuit 
complexity.  When  other  practical  considerations  are  taken  into  account,  the 
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SLEB  technique  has  some  merit  over  other  electron  beam  systems.  This  method 
relies  on  a  fixed  line  geometry  beam,  below  which  samples  are  translated 
mechanically.  The  requirements  for  scanning  electronics  are  eliminated  as 
are  concerns  for  maintenance  of  precise  focusing.  On  an  industrial  scale, 
wafer  processing  speed  should  be  no  slower  than  by  pulsed  beam  processes, 
since  mechanical  translation  will  ultimately  limit  any  automated  beam  process 
Waiting  time  for  system  charging  is  also  eliminated  since  the  beam  will  run 
in  continuous  fashion. 

In  this  work,  I  have  explored  in  some  detail  the  ability  of  the 
SLEB  technique  to  reduce  residual  defects  in  amorphous  layers  and  to  activate 
ion-implanted  dopants  and  traps.  It  will  be  demonstrated  that  future  in¬ 
dustrial  requirements  for  high  quality  junctions  are  indeed  met  by  this 
technique. 

I  do  not  contend  that  pulsed  melt  phase  processes  will  be  without 
future  application.  E-beam  pulsed  diffusion  or  implanted  layer  annealing 
may  prove  important  for  majority  carrier  devices  where  defects  are  not  of 
concern.  Since  photon  energy  absorption  can  be  extremely  shallow,  the  entire 
area  of  contacts  remains  fertile  ground  for  PL  techniques.  However,  at  the 
relative  power  densities  and  speeds  at  which  lasers  can  be  directed,  I  feel 
it  unlikely  that  rastered  laser  techniques  will  ultimately  be  used  to  process 
devices  of  substantial  area.  Even  scanned  CW  laser  systems  require  as  much 
as  10  minutes  to  anneal  a  3  inch  wafer  [61].  Even  longer  cycle  times  should 
be  expected  of  rastered  beam  systems  which  have  significant  charging  times 
between  pulses. 


2.  EXPERIMENTAL  PROCEDURES 

In  this  section,  the  apparatus  and  procedures  used  for  annealing, 
implanting  and  characterizing  semiconductor  material  will  be  discussed  in 
detail. 

2.1  Swept  Line  Electron  Beam  System 

The  annealing  apparatus  is  shown  schematically  in  Figure  2.1.  The 
electron  gun  used  was  an  Alloyed  General  Vacuum  Model  EG3-101  with  modified 
strip  filament.  This  gun  produced  a  stationary,  cyllndrlcally  symmetric 
electron  beam.  An  electrically  grounded  molybdenum  slit  was  positioned  2.5 
cm  below  the  base  of  the  gun  magnet.  The  beam  was  allowed  to  strike  the 
slit  in  an  unfocused  condition.  This  arrangement  produced  a  nearly  rectangular 
electron  beam  spot  on  the  grounded  target  plate.  The  shape  of  the  beam  spot 
could  be  adjusted  through  a  combination  of  gun  magnet  setting,  slit  width  and 
working  distance  to  the  target  plate.  Normally,  a  spot  size  of  1.0  x  1.5  cm 
was  used. 

Samples  were  affixed  to  the  graphite  coated  molybdenum  target  plate 
with  graphite  adhesive  suspension.  This  plate  was  attached  to  a  9  kg  copper 
heat  sink  and  placed  on  a  motor  driven  X-Y  translation  table.  Samples  to  be 
annealed  under  a  common  set  of  conditions  were  arranged  in  a  row,  with  the 
long  dimension  parallel  to  the  long  dimension  of  the  beam  spot.  Each  row  of 
samples  was  annealed  in  a  single  continuous  motion  of  the  X-Y  table.  The 
entire  arrangement  was  enclosed  in  a  stainless  steel  vacuum  chamber.  System 
pressure  was  reduced  to  5  x  10  ^  torr  prior  to  annealing. 

The  electron  beam  was  monitored  by  a  trio  of  Faraday  cups  shown  in 


Figure  2.2.  These  were  arranged  in  a  triangular  array  on  the  X-Y  table.  One 
side  of  this  triangle  was  parallel  to  the  longitudinal  axis  of  the  beam 
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Schematic  diagram  of  the  Swept-Line  Electron  Beam  Annealing  apparatus  (not  to  scale). 
Slit  and  target  plate  are  at  electrical  ground.  Alignment  samples  are  used  to  center 
each  row  of  experimental  samples  on  the  electron  beam.  (See  Yu  |9l|.) 
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spot.  During  annealing,  the  array  was  translated  through  the  beam  spot. 
Collected  cup  currents  were  monitored  continuously  with  a  Keithley  Model 
410CR  electrometer  and  strip  chart  recorder.  These  records  were  used  to 
obtain  beam  energy  density  profiles  and  table  translation  speed. 

Each  cup  was  designed  to  capture  nearly  all  incident  electrons  [62]. 
Secondary  and  backscattered  electrons  generated  within  the  cups  could  have 
have  introduced  considerable  measurement  error.  The  secondary  electrons, 
normally  defined  as  those  with  energy  less  than  50  eV,  [52,63]  were  easily 
captured  by  a  small  attractive  potential  applied  to  the  cup  cores.  The  back- 
scattered  electrons,  however,  could  have  possessed  a  significant  fraction  of 
the  incident  beam  energy  [64].  With  rough  coated  graphite  for  core  material, 
backscattered  electron  yield  was  reduced  to  less  than  8%  of  the  incident 
flux  [52].  The  re-entrant  geometry  of  the  cores  ensured  capture  of  most  of 
those  actually  generated.  I  estimate  no  more  than  2%  of  the  incident  elec¬ 
trons  escaped  capture. 

2.2  Ion  Implantation 

All  implantation  was  performed  on  the  CSL  Accelerators  Inc.  model 
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300-MP,  shown  schematically  in  Figure  2.3.  SI  and  B  F  ions  were 

generated  by  cold  cathode  discharge  from  SiF^  and  BF^  gas  sources  respectively 

Previous  studies  have  shown  that  the  Si  beam  created  with  this  source  will 

contain  less  than  10%  residual  [65].  This  species  cannot  be  removed  by 

mass  analysis.  A  considerably  larger  percentage  of  Nj  could  be  expected  if 

9  + 

the  hot  cathode  and  SiH^  gas  source  were  used.  Be  ions  were  also  generated 

with  the  cold  cathode,  their  source  being  a  solid  Be  canal  insert  which  is 
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continuously  sputtered  by  BF^  gas  discharge.  N  ions  for  GaAsP  implanta¬ 
tion  were  generated  by  hot  cathode  discharge  of  Nj  gas. 


Schematic  diagram  of  the  Accelerators  Inc.  MP-300  ion-implantation  system.  (After  Rosa|l22|.) 


All  samples  were  mounted  7°  off  beam  normal  to  reduce  channeling 
effects  during  implantation.  GaAsP  samples  were  mounted  in  the  target  chamber 
with  conductive  paint.  Silicon  samples  were  mounted  by  direct  pressure 
contacts  to  avoid  contamination  associated  with  the  adhesive.  Si  samples 
amorphized  by  self-implantation  were  mounted  on  a  liquid  nitrogen  cooled 
finger.  The  front  surface  temperature  of  these  samples  was  monitored  with 
a  thermocouple  probe. 

Perfluorinated  polyether  diffusion  pump  fluid  was  used  in  the  target 
chamber  during  all  implantations.  This  fluid  was  found  to  reduce  surface 
carbon  contamination  of  implanted  material  [65].  This  was  particularly  im¬ 
portant  for  low  temperature  implantation,  in  which  residual  gases  condense 
on  the  cold  sample. 

2.3  Methods  of  Characterization 

2.3.1  SEM  channeling  pattern  analysis 

Determination  of  lattice  orientation  and  qualitative  evaluation 
of  crystal  perfection  were  made  with  scanning  electron  channeling  pattern 
analysis.  A  JOEL  JSM-35C  scanning  electron  microscope  was  used.  This  system 
was  equipped  with  Selective  Area  Channeling  Pattern  (SACP)  electronics. 
Electrons  normally  used  for  imaging  were  analyzed  for  Bragg  backscattering. 

The  resulting  patterns  represented  crystal  characteristics  only  for  the  area 
analyzed  by  the  electron  spot.  Since  this  spot  could  be  as  small  as  10  urn 
diameter,  tiny  crystallites  and  disordered  regions  could  be  detected. 

When  operated  in  the  SACP  mode,  the  microscope's  electron  beam  is 


altered  to  produce  a  highly  parallel  beam  [66]  (See  Figure  2.4a).  The  beam 
is  rocked  and  swept  about  a  fixed  point  as  shown  in  Figure  2.4b.  If  the 
specimen  is  crystalline,  the  condition  for  Bragg  backscatter  will  be  satisfied 


Fig  2.4.  (a)  Schematic  diagram  of  the  JEOL  JSM  35C  SEM  beam 

optics  with  channeling  pattern  electronics  activated. 
Field  limiting  aperture  and  lower  deflection  coils  are 
added,  (b)  The  resulting  beam  rocking  is  shown  here 
highly  exaggerated,  (c)  Variation  in  the  first  Bloch 
wave  backscattered  electron  current.  A  large  change 
in  this  current  occurs  when  the  Bragg  condition  is 
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at  certain  points  in  each  scan  described  by: 

2d  sin  9b  »  n  Ae  n  ■  1,2,3,....  (2.1) 

Here,  d  is  the  crystal  plane  spacing,  Aft  is  the  electron  wavelength  and 
9^  is  the  Bragg  angle  (see  Figure  2.4c).  The  resulting  variation  in  back- 

scattered  electron  signal  is  displayed  in  conventional  video  fashion. 

\ 

The  crystallographic  orientation  of  each  sample  area  can  be  deter¬ 
mined  from  its  characteristic  channeling  pattern  [67].  Also,  the  spacings 
of  lattice  planes  may  be  determined  from  the  width  of  the  interference  pat¬ 
tern  lines  as  given  by: 

d  -  cA  /D  (2.2) 

e 

where  D  is  the  SACP  line  width  and  c  is  a  constant  for  fixed  SACP  mode 
magnification  and  electron  energy  [66]. 

One  would  expect  the  quality  of  a  SACP  to  indicate  the  degree  of 
a  sample's  crystal  perfection.  This  is  clearly  shown  in  Figure  2.5,  which 
is  a  succession  of  SACPs  recorded  at  different  points  across  an  amorphous- 
crystal  interface  created  by  partial  annealing.  As  the  analyzing  spot  is 
moved  across,  the  interface  (Figure  2.5a  to  2.5d),  the  channeling  pattern 
lines  related  to  the  low  order  crystal  planes  become  visible.  In  the  figure, 
these  are  the  lines  which  appear  most  nearly  vertical.  Finer  detail  appears 
as  the  beam  moves  onto  undamaged  material.  The  compositional  sensing  mode 
of  the  SEM  is  used  here,  which  tends  to  wash  out  the  horizontal  interference 
lines. 

It  has  been  shown  that  an  amorphous  film  as  thin  as  200  %,  covering 
an  otherwise  crystalline  sample,  will  introduce  sufficient  Interference  to 
wash  out  a  SCAP  [67].  Physical  surface  defects  such  as  scratches,  stress 


2. 3. 2.1  Overview  and  apparatus  description 


Deep  Level  Transient  Spectroscopy  (DLTS)  was  used  in  this  work  to 
detect  and  analyze  semiconductor  defects.  The  theoretical  basis  for  DLTS 
has  been  well  established:  [68,69]  however,  I  will  briefly  review  the  rele¬ 
vant  fundamental  concepts.  Basically,  DLTS  operates  through  analysis  of 
transient  capacitance  decay  of  a  reverse  biased  diode.  An  example  for 
majority  carrier  traps  in  the  depletion  layer  of  a  p+  -  n  diode  is  illustrated 
in  Figure  2.6.  The  device  under  test  is  held  at  a  quiescent  reverse  bias. 

The  traps  are  initially  free  of  electrons.  A  short  bias  reducing  pulse  is 
applied  which  compresses  the  depletion  region,  allowing  the  traps  to  fill. 
Immediately  after  the  pulse  ends,  the  junction  capacitance  falls  beloy  its 
quiescent  value  C,  by  an  amount  AC  due  to  the  compensating  charge  of  the 
filled  traps.  Electrons  are  then  thermally  emitted  from  the  traps,  and  the 
capacitance  returns  to  its  quiescent  value.  The  characteristic  time  of  the 
capacitance  decay  depends  primarily  upon  the  ratio  of  the  trap  activation 
energy  and  device  temperature.  By  measuring  this  decay  time  as  a  function  of 
temperature,  the  energy  level  of  the  defect  may  be  inferred.  The  magnitude 
of  the  capacitance  change  AC  is  a  measure  of  the  concentration  of  defects. 

By  adjusting  the  magnitude  of  the  bias  reducing  pulse,  one  may  determine  the 
distribution  of  a  particular  defect.  The  trap  capture  cross  section  can  be 
obtained  from  a  measurement  of  AC  with  varying  bias  reducing  pulse  width. 

If  the  bias  reducing  pulse  magnitude  is  increased  so  as  to  cause  actual 


injection,  DLTS  will  detect  minority  carrier  traps  in  the  same  way.  The 
sign  of  AC  determines  whether  majority  or  minority  carrier  traps  are  being 
analyzed.  It  can  be  seen  that  DLTS  is  a  powerful  technique  for  the  detection 
and  analysis  of  semiconductor  defects. 

The  DLTS  system  used  in  this  work  was  developed  by  D.  S.  Day  et  al. 
[70]  and  is  shown  schematically  in  Figure  2.7.  This  system  employs  a  unique 
two-arm  bridge  circuit.  Each  arm  contains  a  diode  of  similar  C-V  and  I-V 
characteristic.  These  are  mounted  on  a  common  TO-18  header  and  held  within 
a  variable  temperature  dewar.  This  arrangement  keeps  differences  in  device 
temperature  to  a  minimum.  The  bridge  is  driven  by  a  20  MHz  modulating  signal. 
The  device  labeled  "Test  Diode"  is  driven  exactly  180°  out  of  phase  with  the 
"Dummy  Diode".  Shifts  in  the  characteristics  of  the  test  device  tend  to  be 
cancelled  by  the  dummy  device,  and  the  entire  bridge  remains  in  balance 
during  temperature  scanning.  Note  that  the  DC  quiescent  bias  is  applied  to 
both  devices,  while  only  the  test  device  receives  the  bias  reducing  pulse. 
Isolation  from  this  transient  pulse  is  provided  by  two  attenuators,  a  power 
splitter,  a  phase  shifter  and  bandpass  filter.  We  may  therefore  assume  that 
only  the  test  device  produces  a  AC  signal. 

The  bridge  signal  is  demodulated  at  the  HP  mixer,  amplified  and 
finally  detected  at  the  lock-in  amplifier.  A  special  gate  circuit  eliminates 
the  large  capacitance  overshoot  present  during  the  bias  reducing  pulse. 

This  prevents  amplifier  overload  and  keeps  system  response  time  limited  to 
less  than  1  usee.  The  magnitude  of  AC  is  displayed  on  the  X-Y  recorder  as  a 
function  of  device  temperature.  A  Booton  76-3A  standard  capacitor  is  used  in 
place  of  the  diodes  to  align  the  bridge  to  the  capacitive  mode  and  to  provide 


calibration. 
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Fig.  2.7.  Block  diagram  of 

the  two- 

-diode  DLTS  System. 

(After  Day|70|) 

2. 3. 2. 2  DLTS  data  analysis 

In  this  work,  only  the  energy  level  of  the  traps  and  their  concen¬ 
tration  profiles  were  studied  in  detail.  The  manner  in  which  this  informa¬ 
tion  was  obtained  from  DLTS  data  is  summarized  in  this  section  and  in  Appen¬ 


dix  A. 


From  the  principle  of  detailed  balance  and  carrier  statistics,  the 


emission  rate  of  carriers  from  a  trap  (e),  can  be  expressed  in  terms  of  its 
capture  cross  section  (a),  the  trap  degeneracy  factor  (g),  the  rms  thermal 
velocity  of  the  captured  carrier  (<v>),  the  appropriate  effective  density 
of  states  (N)  and  the  depth  of  the  trap  in  the  bandgap  (AE)[69].  For  the 
case  of  electron  traps  (and  adding  appropriate  subscripts) ; 


a  <v  >N 
n  n  c 


exp  (-AE/kgT) 


(2.3) 


If  one  assumes  no  temperature  dependence  of  the  capture  cross  section. 


(valid  for  certain  types  of  traps)  and  using: 


N  -  2M 
c  c 


f2irm  kgT  |  3/2 


(2.4a) 
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<v  >  * 
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3kgT  )  1/2 


(2.4b) 


we  may  write; 


±  -  e  «  T2.  exp CAEV) 
n 


(2.5) 


where  x^  is  the  time  constant  corresponding  to  this  emission  rate.  Finally, 


T2xn  oc  exp  (AE/kBT) 


(2.6) 


Thus  a  plot  of  ln(T  x^)  vs  1/T  will  yield  a  line  whose  slope  is  proportional 
through  known  constants  to  AE.  x  is  determined  by  the  fundamental  frequency 


set  on  the  lock-in  amplifier.  T  is  the  temperature  at  which  a  trap  pro¬ 
duces  the  largest  AC  signal.  Each  trap  can  be  analyzed  at  several  different 
frequencies,  and  least  squares  fit  to  these  data  can  then  be  made  to  improve 
accuracy  [71]. 

The  actual  magnitude  of  the  capacitance  transients  AC  must  be  cal¬ 
culated  from  the  output  of  the  lock-in  amplifier.  This  output  is  the  time 
average  of  the  product  of  the  input  signal  and  the  lock-in's  square  wave 
weighting  function.  Both  the  shape  of  the  input  waveform  and  its  phase  with 
respect  to  the  weighting  function  will  affect  the  value  of  AC  detected.  In 
this  work  we  follow  the  analysis  of  Day  [71,72],  The  input  waveform  is 
assumed  to  be  a  train  of  simple  exponentially  decaying  spikes.  The  weighting 
function  is  adjusted  to  be  in  phase  with  the  leading  edge  of  the  bias  re¬ 
ducing  pulse  (bias-pulse  phase  reference  method) .  Under  this  condition, 
mathematical  analysis  of  the  detection  process  yields  the  expression: 

AC  ■  0.560  x  x  Sg  x  (2.7) 


where  AC  is  the  transient  capacitance  in  pF,  SL  is  the  lock-in  sensitivity 
in  volts,  Sn  is  the  calibration  factor  in  pF/V  and  V  is  the  lock-in  output 
in  volts. 

From  these  values  of  AC,  one  may  construct  a  profile  of  defect 
density.  If  one  considers  the  example  of  a  p+-n  diode,  Lang  has  determined 


*  #  * 
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(2.8) 


N  (x  ) 
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where  C  and  W  correspond  to  the  diode  capacitance  and  depletion  width  while 
quiescent  reverse  bias  V  is  applied.  5  O^r)  is  the  change  in  the  ~  ratio 
between  two  successive  measurements,  with  bias  reducing  pulse  and 


V  +  6V  .  N  is  the  fixed  charge  concentration  (due  to  ionized  donors  in 
c  c 

our  example)  at  the  n  edge  of  the  depletion  layer.  N„(x  )  and  N+(x  )  are 

T  C  C 

respectively  the  densities  of  traps  and  of  fixed  charge  at  the  location 
(xc)  of  the  depletion  layer  edge  during  the  bias  reducing  pulse  (v^  +  AV^) . 

A  trap  profile  may  be  constructed  via  a  series  of  measurements  of  AC  while 
varying  Vc>  with  fixed  V.  To  be  valid,  AC«C  and  traps  must  occur  only  on 
the  n-side  of  the  junction.  It  must  also  be  assumed  that  the  acceptor  and 
donor  concentrations  change  in  the  same  ratio  on  both  sides  of  the  metal¬ 
lurgical  junction.  This  assumption  is  valid  for  linearly  graded  and  step 
junctions. 

The  fixed  charge  profiles  were  produced  via  point-by-point  capaci¬ 
tance  voltage  measurements.  A  Hewlett  Packard  6114A  precision  power  supply 
provided  bias  to  a  Boo ton  72B  capacitance  meter  for  these  measurements.  A 
sufficiently  small  interval  in  voltage  was  placed  between  measurements  to 
estimate  dC/^.  The  magnitude  of  fixed  charge  concentration  was  determined 
via  the  expression: 

3 

N(x) - -  (2.9) 

eqA  (-dC(V) / dV) 

The  depth  x  was  determined  by  numerical  integration  of  the  fixed  charge 
concentration  profile  and  knowledge  of  the  device  geometry  (See  Appendix  A, 
Section  A. 2,  lines  19800-21800). 

Trap  profile  analysis,  including  determination  of  the  fixed  charge 
profile,  was  accomplished  by  computer.  The  FORTRAN  program  TRAPSI  was 
developed  specifically  for  Schottky  barrier  diodes  on  n-type  silicon.  This 
code  is  described  in  detail  in  the  appendix. 


2*3.3  Differential  resistivity  and  Hall  effect  measurements 


Profiles  of  electrically  active  implanted  dopants  were  established 
by  differential  resistivity  and  Hall  effect  measurements  made  in  conjunction 
with  chemical  layer  stripping.  The  apparatus  was  developed  by  McLevige 
et  al.  [74],  and  is  shown  schematically  in  Figure  2.8.  This  double  AC  Hall 
system  possessed  high  sensitivity,  good  noise  rejection  and  inherent  averaging 
of  misalignment  and  thermoelectric  effects.  During  Hall  measurements,  the 
samples  were  placed  in  an  AC  electromagnet  which  was  driven  at  one  fourth  the 
basic  system  frequency  (f^ )  of  1  kHz.  Via  a  mixer  and  filter,  a  s*-£nal 

was  generated  which  formed  the  lock-in  reference  signal  for  the  Hall  voltage, 
f^  was  used  as  the  lock-in  reference  for  resistivity  measurements. 

Samples  were  prepared  by  the  technique  described  by  Tsai  [65].  The 
carrier  type  of  the  substrate  was  choosen  so  that  a  p-n  junction  was  formed 
after  implantation  and  annealing.  First  a  120  mil  square  van  der  Pauw  mesa 
pattern  was  photolithographically  defined  and  plasma  etched.  Gold  contacts 
were  then  evaporated  and  later  sintered  at  300°C  in  flowing  dry  H^,.  Samples 
showing  good  junction  characteristics  were  selected  so  that  only  the  im¬ 
planted  layer  contributed  to  the  measurement.  After  each  measurement,  a 
thin  layer  of  the  sample  was  removed  by  chemical  etching.  A  series  of  wax 
defined  mesas  were  formed  on  each  sample  during  thinning.  The  heights  of 
these  mesas  were  later  measured  by  Detak  mechanical  stylus  to  establish  the 
depth  scale.  The  thinning  and  measurement  process  was  continued  until  the 
samples  became  too  resistive  to  measure. 

The  sheet  resistivity  (ps),  and  sheet  Hall  coefficient  (Rg) ,  are 


calculated  from  [75]: 
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Fig  2.8.  Block  diagram  o£  the  double  AC  Hall  Apparatus. 
(After  McLevlge  et  al.  I  74 1  ) 
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s  ^  tVBDAC<m8)  • 
3  B(rms)  I(rms) 


(2.9) 


(2.10) 


Here  the  designation  corresponds  to  the  voltage  between  contacts  A  and 

B  with  current  I  passing  between  contacts  C  and  D,  and  similarly  for  V___A. 

BCDA 

A  VBDAC  is  the  chan8e  in  Hall  voltage  induced  by  placing  the  sample  into  the 

magnetic  field  of  strength  B.  The  correction  factor  f  (  ■= -  )  has  been 

BCAD 

tabulated  by  van  der  Pauw  [75]. 

The  average  mobility  and  carrier  concentration  of  a  particular 
layer  j  of  thickness  d^  are  given  by:  1 76 1 


n  A  A  R 

u  -  "I  t  -  81-1  81-1  ~  8i 
HJ  p8j  p8J-l  p8j-l  "  p8j 


(2.11) 
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8J  **J-1 

edjWi 


(2.12) 


The  designation  j-1  refers  to  the  preceedlng  layer,  y  is  the  ratio  of  Hall 

(yu  )  and  conductivity  mobility  (u.).  Its  value  was  taken  to  be  0.73,  which 
Hj  J 

is  appropriate  for  the  highly  doped  p-type  silicon  layers  studied  here  [77]. 


3.  SLE2  ANNEALING  OF  Si  -IMPLANTED  AMORPHOUS  SILICON 

This  chapter  presents  a  comparative  study  of  amorphous  ion-implanted 
silicon  recrystallized  by  svept-line  electron  beam  and  conventional  furnace 
annealing.  Solid  phase  recrystallization  of  amorphous  silicon  has  been 
studied  previously  using  furnace  techniques  j 65, 78, 79 |.  Dopants  implanted 
in  high  doses  can  be  efficiently  activated  with  accompanying  recrystalliza- 
tlon  at  relatively  low  temperatures  (approximately  550°C).  As  demonstrated 
by  Tsai  et  al.  | 80 { ,  however,  electrical  profiles  of  these  dopants  can  in¬ 
clude  an  electrically  in  active  tall  due  to  heavy  residual  implantation 
damage.  (See  section  4. 1)  This  layer  is  located  between  the  bulk  material 
not  damaged  by  the  original  implantation  and  material  which  regrows  completely 
during  annealing.  In  a  p-n  junction,  this  damage  coincides  with  the  de¬ 
pletion  region,  and  will  adversely  affect  leakage  currents  and  minority  car¬ 
rier  lifetimes. 

Tsai  et  al  1 81]  have  also  shown  that  the  location  of  the  residual 
damage  can  be  pushed  well  beyond  the  p-n  junction  by  amorphizing  the  sample 
(eg.,  by  Si+  ion  implantation)  to  a  depth  beyond  the  implanted  dopant  profile. 
Samples  prepared  in  this  manner  have  completely  active  dopant  profiles.  In 
fact,  fluorine  migration  to  defects  | 81 |  has  shown  that  the  unannealed  damage 
is  moved  to  the  tail  of  the  deep  amorphizing  implantation. 

The  presence  of  this  damage  was  further  confirmed  through  a  pre¬ 
liminary  study  using  silicon  solar  cells.  This  work  was  undertaken  in 
cooperation  with  Dr.  J.  Eguren  of  the  Instituto  De  Energia  Solar,  Madrid, 

Spain.  A  set  of  n+pp+  back  surface  field  cells  was  fabricated  with  and 
without  the  use  of  amorphizing  self-implantation.  Here  the  low-hi  (pp+) 
region  is  °eant  to  , reduce  the  back  surface  recombination  velocity.  The 


front  surface  n+  layers  were  formed  by  conventional  furnace  diffusion. 

Some  of  the  samples  then  underwent  multiple  low  temperature  Si  implantation 
to  produce  a  continuous  amorphous  back  surface  layer.  A  single  room  tempera¬ 
ture  gallium  implant  was  then  made  in  all  the  samples  to  form  the  back  sur¬ 
face  p+  layer.  The  energy  of  this  implantation  was  sufficiently  low  that 
all  gallium  atoms  remained  well  within  the  amorphous  layer.  After  furnace 
annealing,  the  solar  cells  were  characterized  under  Air  Mass  1  Illumination. 

Representative  results  are  shown  in  Table  3.1.  These  data  show 

consistently  better  performance  in  cells  which  were  not  pre-amorphized. 

Samples  1.2  -  1.4  show  generally  smaller  open  circuit  voltage  (V  )  and 

ocf,b 

short  circuit  currents  (Ij  when  compared  with  sample  1.1.  Implantation 

and  annealing  conditions  were  such  that  one  would  expect  a  completely  active 
gallium  layer  in  the  pre-amorphized  cells  and  only  partial  activation  of 
gallium  in  the  unamorphized  cells.  This  alone  would  lead  one  to  expect 
just  the  opposite  of  the  observed  difference  in  performance.  The  result  can 
be  explained  in  terms  of  reduced  cell  base  lifetime,  due  to  the  self-implanta¬ 
tion  induced  damage.  The  study  described  below  will  characterize  this  damage 
in  some  detail. 

DLTS  is  used  here  to  directly  evaluate  residual  crystal  damage. 

The  experiment  is  simplified  through  the  use  of  Schottky  barrier  diodes 
fabricated  on  recrystallized  self-implanted  material.  Ambiguities  of  p-n 
junction  dopant-defect  complexes  are  thus  eliminated.  Since  Schottky  diodes 
can  be  prepared  without  disturbing  the  bulk  crystal,  a  direct  comparison  of 
SLEB  and  furnace  recrystallized  material  can  be  made. 


3 . 1  Experimental  Procedure 


Devices  used  in  this  study  were  fabricated  from  0.24  ft-cm 
16  -3 

(4  x  10  cm  ) ,  phosphorus  doped  Czochralski  grown  bulk  silicon.  This  level 
of  dopant  was  chosen  to  provide  the  best  compromise  between  Debye  length 
limits  on  profile  resolution  and  the  requirement  for  full  depletion  of  the 
regrown  amorphous  layer  during  measurement  [83].  As  in  the  solar  cell  ex¬ 
periments  described  above,  the  samples  received  four  cold  (<,  -100°C)  amor- 
phizing  implants  with  Si*  ions  of  energies  32,  91,  158  and  255  keV.  These 

energies  and  the  corresponding  doses  were  chosen  to  give  the  best  fit  to 

20  +  3 

an  ideal  0.5  ym  deep,  flat  implanted  profile  of  1  x  10  Si  /cm  .  Etching 
experiments  [81]  determined  the  depth  of  the  amorphous  layer  to  be 
4550  X  -  125  X. 

Samples  were  then  recrystallized  either  by  the  SLEB  process  or  by 
conventional  furnace  annealing.  The  electron  beam  annealing  process  was 
described  in  section  2.1.  An  electron  energy  of  20  keV  and  table  translation 
speed  of  0.32  cm/sec  were  used  in  this  work.  A  second  set  of  samples  was 
thermally  annealed  in  a  quartz-lined  furnace  tube  under  a  forming  gas  ambient. 
Unimplanted  material  was  heat  treated  simultaneously  to  detect  defects  which 
might  be  Introduced  during  furnace  annealing  and  unrelated  to  the  implanta¬ 
tion. 

After  annealing,  the  front  surfaces  of  the  samples  were  anodically 
oxidized  [65]  to  remove  a  200  X  layer  of  Si,  exposing  clean  surfaces  for 
Schottky  barrier  formation.  A1  Schottky  barriers  and  Au-Sb  ohmic  back-side 
contacts  were  formed  by  vacuum  evaporation  through  shadow  masks.  Individual 
devices  were  cut,  mounted  on  headers,  and  thermocompression  lead  bonded. 

Sample  heating  did  not  exceed  300° C  for  10  minutes  during  this  processing. 


Defect  energy  levels  and  profiles,  as  well  as  profiles  of  fixed 

charge  were  measured  and  calculated  by  the  methods  described  in  section  2.3.2 

and  the  appendix.With  the  exceptions  noted  below,  measured  values  of  AC/C 
-2  -5 

were  between  10  and  10  ,  satisfying  the  restriction  that  AC/C«1  168,69]. 

3.2  Results  and  Discussion 

DLTS  spectra  are  presented  in  Figure  3.1  as  a  function  of  30-minute 

isochronal  furnace  annealing.  The  DLTS  signal  is  plotted  logarithmically, 

with  majority  carrier  signal  increasing  along  the  vertical  scale.  The  peak 

trap  signal  decays  as  the  regrowth  temperature  is  increased,  as  expected. 

In  each  case,  the  spectra  are  dominated  by  two  defect  signals.  The  lower 

temperature  signal  (corresponding  to  Ec  -  0.22  eV)  appears  in  all  three 

cases  and  becomes  dominant  as  annealing  temperature  is  increased.  The 

higher  temperature  trap  signal  changes  from  E  -  0.50  eV  to  E  -  0.35  eV  to 

c  c 

E  -  0.42  eV  as  the  annealing  temperature  is  increased.  Defects  with 
activation  energies  within  +  0.02  eV  of  these  four  trap  signals  have  been 
observed  in  implanted  [84],  irradiated  [85,86],  or  laser  damaged  silicon 
[87].  Defect  spectra  for  control  Schottky  barriers  fabricated  on  unimplanted 
but  simultaneously  annealed  Si  show  no  detectable  defect  signal. 

Similar  spectra  are  shown  in  Figure  3.2  for  SLEB  samples  processed 
at  three  peak  beam  power  densities.  There  is  a  marked  reduction  in  the 
peak  trap  signal  from  lowest  to  highest  power  density.  The  dominant  signal 
peak  of  the  well-annealed  SLEB  case  (38  W/cm  )  is  about  a  factor  of  four 
smaller  than  the  smallest  peak  in  the  800°C  case  of  Figure  3.1.  The  designa¬ 
tion  ”x2"  refers  to  a  sample  which  was  annealed  by  two  passes  of  the  beam. 

Notice  that  the  double  annealed  sample  at  Intermediate  power  density  (27.5 
2 

W/cm  )  shows  the  largest  over-all  defect  signal  spectra.  The  thermal 


DLTS  Spectra 


DLTS  spectra  for  Schottky  barriers  fabricated  on  furnace  annealed  silicon.  Majority 
carrier  signal  is  plotted  on  the  logarithmic  vertical  scale.  (See  Soda  1 31 j.) 


spectra  for  samples  annealed  by  SLEB  for  three  power  densities.  The  vertical 
:  has  units  Identical  to  figure  3.1.  (See  Soda  1 31 1 . > 


activation  energies  observed  in  the  SLEB  case,  with  the  exception  of 

E^  -  0.34  eV,  do  not  correspond  to  those  observed  in  the  furnace  case. 

Activation  energies  within  +  0.02  eV  of  the  observed  levels  E  -  0.17  eV, 

—  c 

0.18  eV,  0.34  eV,  and  0.46  eV,  have  been  previously  observed  in  implanted 
[84],  irradiated  [85,86],  or  laser  damaged  [87]  silicon.  The  level  E£  -  0.80 
eV  has  not  been  reported  previously. 

Profiles  of  fixed  charge  and  defect  concentration  for  the  furnace 
annealed  devices  are  presented  in  Figures  3.3  and  3.4.  In  each  case,  the 
profiles  are  dominated  by  a  layer  of  defects  and  fixed  charge  near  the 
original  amorphous-crystalline  interface  (marked  by  the  arrow) .  The  magni¬ 
tude  of  these  effects  decays  with  increased  annealing  temperature.  Also 
observed  is  a  relatively  small  concentration  of  defects  in  the  near  surface 
regions.  This  is  accompanied  by  a  tendency  for  fixed  charge  concentrations 
to  return  to  values  corresponding  to  the  original  bulk  dopant  concentration 
(4  x  10*^  cm  "*) . 

In  the  700°C  and  800°C  cases,  the  defect  profiles  reach  a  maximum 
at  0.42  to  0.45  pm  respectively,  corresponding  to  the  measured  thickness  of 
the  original  amorphous  layer.  This  observation  verifies  the  location  of 
damage  reported  by  Tsai  et  al.  [81,82]  from  their  Hall  profiling  and  fluorine 
decoration  studies.  Previous  studies  in  irradiated  material  also  indicate 
that  the  near  surface  recrystallized  regions  tend  to  have  fewer  defects  [88], 
as  observed  here. 

The  fixed  charge  profiles  presented  here  must  include  the  contribu¬ 
tions  of  both  ionized  donors  and  defects  [83].  In  regions  where  the  defect 
concentrations  become  comparable  to  the  fixed  charge  profiles,  one  cannot 
strictly  interpret  the  latter  as  a  donor  profile.  Rather,  the  fixed  charge 
profile  represents  N.+nNT,  where  n  and  its  sign  depend  upon  the  charge  state 
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Fig.  3.3. 


Fixed  charge  and  defect  profiles  for  700  C  furnace  processed 
silicon.  The  arrow  Indicates  the  location  of  the  original 
amorphous-crystalline  transition  region.  (After  Soda  1311.) 
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of  Che  defect.  This  is  important  in  the  700°C  case  near  the  peak  damage 
region.  In  the  800° C  case,  however,  ionized  donor  charge  predominates  and 
we  see  definate  donor  redistribution,  particularly  in  the  neighborhood  of 
the  original  amorphous-crystalline  interface.  The  donor  redistribution 
effect  can  be  understood  in  terms  of  a  dissociative  diffusion  process  of  the 
type  proposed  by  Baruch  [89].  Motion  of  dopants  is  explained  by  a  pair  of 
competitive  processes  in  which  a  defect  will  first  liberate  a  donor  atom 
from  its  lattice  site,  then  cause  the  donor  to  move  by  attraction  to  another 
defect  site.  The  donor  atom  undergoes  a  net  motion  in  the  direction  opposite 
the  motion  of  defects.  These  processes  are  used  to  explain  bulk  dopant  re¬ 
distribution  and  peaking  in  proton  irradiated  silicon  [89].  The  shapes  of 
the  donor  profiles  in  the  near-surface  region  suggest  dopant  motion  toward 
the  interface  region. 

It  is  likely  that  the  solid  phase  regrowth  process  itself  is  in¬ 
volved  in  the  redistribution  of  both  defects  and  donor  atoms,  especially 
in  the  600° C  case  (not  shown) .  A  reasonable  extrapolation  of  previously 
determined  solid-phase  epitaxial  regrowth  rates  [78]  indicate  that  nearly 
25%  of  the  total  annealing  time  was  required  to  recrystallize  the  amorphous 
layer  in  this  case.  The  advancing  crystal  front  seems  to  push  defects 
ahead  of  it,  accelerating  the  motion  of  dopant  atoms  and  leaving  a  bulge 
in  the  defect  profiles  near  0.23  urn.  Regrowth  times  are  estimated  to  be 
factors  of  14  and  370  times  faster  respectively  for  the  700°  and  800°  anneals 
shown  here.  No  bulge  in  the  defect  profile  is  observed  in  the  regrown 
regions  of  Figures  3.3  and  3.4,  indicating  that  either  the  crystal  front  has 
moved  too  rapidly  to  push  the  defects,  or  they  have  annealed  extensively 
after  regrowth  is  complete. 


material  are  shown  in  Figures  3.5  and  3.6.  In  both  cases,  the  residual 

defect  concentrations  are  small  enough  that  one  may  interpret  the  fixed 

charge  profiles  as  due  entirely  to  ionized  donors.  After  a  single  sweep 
2 

by  a  19  W/cm  electron  beam,  the  sample  exhibits  a  series  of  broad,  rela¬ 
tively  flat  defect  profiles  (not  shown) .  Surface  channeling  pattern 
analysis  shows  that  annealing  in  this  case  is  insufficient  to  completely 
regrow  the  amorphous  layer. 

Recent  studies  of  nitrogen  implanted  GaAs^  1 91  j  have  suggested 

that  multiple  scanning  with  SLEB  can  improve  crystal  quality  without  defect 

2 

migration.  This  possibility  was  explored  in  the  27.5  W/cm  sample  (Figure 
3.5),  which  was  double  annealed.  In  this  case,  severe  dopant  redistribu¬ 
tion  occurs.  Electron  microscope  channeling  analysis  shows  complete  re¬ 
crystallization  was  not  achieved.  Since  the  0.80  eV  defect  AC/C  signal  is 
so  large  (*  0.1),  its  profile  is  doubtful.  Some  surface  dopant  out-diffusion 
is  observed.  A  large  misalignment  of  the  defect  and  donor  peaks  with  the 
amorphous-crystal  Interface  is  observed,  probably  due  to  a  near  surface 
donor  depleted  region  which  cannot  be  analyzed  by  C(V)  measurements.  We 
believe  the  depth  scale  is  affected,  and  that  the  peaking  of  defects  and 
charge  in  Figure  3.5  actually  occurs  at  the  original  interface  (marked  by 
the  arrow) . 

The  threshold  for  complete  crystal  growth  is  crossed  in  the  38 

2 

W/cm  case  (Figure  3.6).  SEM  analysis  shows  perfect  near  surface  crystal 
order.  Both  dopant  redistribution  and  out-diffusion  effects  are  minimized. 
The  profile  of  the  single  defect  observed  in  this  case  (E£  -  0.32  eV)  occurs 
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at  about  the  depth  of  the  original  amorphous-crystalline  transition.  Peak 
defect  concentrations  are  a  factor  of  ten  smaller  than  in  the  800° C  furnace 
annealed  case,  while  near  surface  defect  densities  are  a  factor  of  five  to 
ten  smaller. 

Other  studies  of  recrystallization  of  amorphous  Si  by  pulsed  E-beam 
processes  [92]  explain  the  regrowth  strictly  in  terms  of  beam  induced  heating. 
In  this  work,  total  annealing  times  on  the  order  of  two  seconds  were  used. 
These  are  considerably  longer  than  either  pulsed  [92]  or  raster  scanned  [93] 
E-beam  processes.  Therefore,  the  annealing  conditions  are  somewhat  more 
similar  to  those  observed  to  cause  MeV  electron  or  keV  proton  irradiation 
induced  defect  annealing.  In  these  studies,  TEM  techniques  [94,95]  and  IR 
absorption  [96]  are  used  to  observe  both  crystal  regrowth  and  accelerated 
defect  annealing. 

Geranismenko  [97]  has  suggested  that  point  defects  generated  by 
higher  energy  electrons  can  be  responsible  for  decay  of  relatively  stable 
defect  formations,  such  as  those  observed  in  this  work.  These  defects  are 
thought  to  compensate  or  induce  migration  of  the  defects  responsible  for 
the  original  amorphous  state.  He  also  suggests  that  the  electron  induced 
ionization  may  change  the  charge  state  of  the  defects,  increasing  their 
mobility. 

It  would  seem  that  generation  of  point  defects  by  electrons  of 
only  20  keV  is  a  minor  effect.  The  lower  limit  for  electron  induced  atomic 
displacements  is  145-200  keV  [98,99].  However,  Hinkley  [100]  has  invoked  a 
two-step  energy  transfer  process  to  explain  subthreshold  defect  generation 
in  silicon.  In  this  process,  incident  low  energy  electrons  excite  valence 
electrons,  producing  plasmons.  These  plasmons  transfer  their  energy  to  the 


lattice,  eventually  causing  atomic  displacements.  The  plasmon  process  also 
explains  an  experimentally  observed  [100]  reduction  in  this  threshold  with 
increasing  sample  temperature.  Since  the  SLEB  produces  sample  heating,  the 
plasmon  process  may  contribute  to  the  observed  reduction  in  residual  defects, 
such  as  observed  in  Figure  3.6. 

Some  support  for  this  theory  comes  from  comparison  of  peak  trap 

concentrations  of  the  beam  annealed  samples.  If  a  strictly  thermal  regrowth 

and  annealing  process  were  responsible,  one  would  observe  a  steady  overall 

reduction  in  the  trap  concentration  as  beam  power  is  increased.  Double 

annealing  should  accelerate  this  trend.  In  fact,  the  largest  peak  trap 

concentrations  and  significant  donor  redistribution  are  observed  in  the 
2 

27.5  W/cm  double  sweep  case,  which  was  certainly  warmer  during  the  second 
beam  sweep.  In  this  case,  plasmon  induced  defects  could  explain  the 
accelerated  donor  migration.  No  major  donor  redistribution  effects  are 
observed  in  samples  treated  with  a  single  sweep. 

3.3  Conclusions 

The  data  show  that  relatively  deeply  amorphlzed  Si  layers  can  be 
recrystallized  by  swept-line  electron  beam  processing.  DLTS  analysis  shows 
that  SLEB  can  reduce  residual  defect  densities  to  levels  well  below  that 
achievable  by  conventional  furnace  annealing.  This  is  particularly  time  in 
the  transition  region  between  bulk  and  recrystallized  material,  where  residual 
damage  is  resistant  to  high  temperature  annealing.  It  is  significant  that 
SLEB  recrystallization  can  be  effective  without  the  use  of  additional  furnace 
treatment.  SLEB  processes  can  also  reduce  bulk  dopant  redistribution  effects 
which  occur  during  furnace  regrowth.  These  results  do  not  represent  the 


4.  SLEB  ANNEALING  0?  BF^  IMPLANTED  SILICON 
This  chapter  extends  the  study  of  SLEB  recrystallized  amorphous 
silicon  to  Include  BF^  implanted  Si,  in  which  the  amorphization  is  created 
by  the  fluorine  and  acceptor  doping  is  provided  by  the  boron.  Hall  and 
Secondary  Ion  Mass  Spectroscopy  (SIMS)  profiling  reveal  high  electrical 
activation  efficiency  and  only  modest  spatial  redistribution.  Unlike  furnace 
annealed  material,  active  boron  is  observed  in  the  as-implanted  amorphous- 
crystalline  transition  region.  This  can  be  important  to  the  production  of 
high  quality  implanted  p-n  junctions.  Distribution  effects  of  Implanted 
fluorine  are  also  explored,  and  these  migration  effects  suggest  that  SLEB 
annealing  does  not  produce  crystal  damage  beyond  the  original  implanted 
profile. 

4.1  Experimental  Procedure 

The  samples  used  in  this  study  were  Czochralski  grown,  phosphorus 

doped,  0.24  fl-cm,  <100>  silicon.  Implantation  of  BF^  ions  was  performed 
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at  150  keV  to  a  fluence  of  1  x  10  cm  at  room  temperature.  These  condi¬ 
tions  have  been  found  sufficient  to  create  an  amorphous  layer  ~  1325  X 
thick  [80].  Beam  annealing  was  accomplished  as  described  in  section  2.1. 
Samples  of  dimensions  1  x  1. 3  cm  were  attached  to  the  sample  stage  and 
translated  at  0.32  cm/sec  during  annealing.  An  electron  energy  of  20  keV 
was  used  in  this  study. 

After  annealing,  gold  contacts  were  applied  and  a  van  der  Pauw 
geometry  mesa  was  defined  by  photolithographic  techniques  and  plasma  etching. 
Each  slice  accomodated  either  six  or  eight  Hall  patterns.  Heat  cycling 
after  SLEB  annealing  did  not  exceed  300°C  for  15  sec.  Electrical  carrier 
profiling  was  accomplished  by  double  AC  Hall  measurements  in  conjunction  with 


chemical  layer  stripping,  as  described  in  section  2.3.3.  Etch  depths  were 
measured  at  various  points  in  the  profile  by  using  a  Sloan  Dektak  mechanical 
stylus. 


Samples  used  in  electrical  profiling  measurements  are  subsequently 
used  in  SIMS  profiling,  to  facilitate  comparisons  between  atomic  and  electri* 
cally  active  profiles.  Gold  contacts  were  removed  exposing  processed  im¬ 
planted  silicon  which  was  uneffected  by  the  electrical  characterization.  A 
maximum  of  65  mils  separated  the  regions  which  were  electrically  and 
chemically  profiled.  Secondary  ion  mass  spectrometry  was  then  performed 
using  the  Cameca  AMS  3F  ion  microprobe.  0*  ions  of  14.5  keV  were  employed 
for  sputtering.  Mass  11,  19,  28,.  and  31  signals  were  monitored  during  each 
profile  by  peak  switching  techniques.  Concentration  scales  were  established 
through  reference  to  an  unannealed  standard.  A  Tencor  mechanical  stylus  was 
used  to  measure  crater  depths  and  thus  establish  the  depth  scale. 

4. 2  Results 

SLEB  annealing  was  performed  at  power  densities  of  between  19  and 

2 

64  W/cm  .  SEM  channeling  pattern  analysis  indicated  that  all  samples 

annealed  in  this  range  attained  high  quality  near  surface  crystal  order. 

Those  annealed  at  the  lowest  power  density,  however,  were  too  resistive  to 

be  Hall  profiled.  The  net  active  carrier  profiles  for  the  three  electron 

beam  power  density  conditions  are  shown  in  Figure  4.1.  The  designation  "X3" 

means  that  the  sample  was  passed  below  a  beam  of  the  indicated  power  density 

three  times.  A  ninety  second  interval  was  allowed  between  successive  sweeps 

under  the  beam.  We  observe  a  general  increase  in  the  peak  active  boron  con- 

2 

centration  as  the  annealing  power  density  Increases  from  30  to  64  w/cm  . 

2 

The  64  w/cm  profile  compares  favorably  with  the  as-implanted  SIMS  atomic 
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Fig  4.1.  Ne£  acceptor  concentration  and  conductivity  mobility  profiles  for 
BF2  implanted  silicon  annealed  at  three  SLEB  power  densities. 
Shown  for  comparison  are  the  as-implanted  atomic  boron  and 
corresponding  "ideal"  mobility  profile.  The  furnace  annealed 
profiles  are  due  to  Tsai  et  al.  | 80 | .  The  calculated  energy 
deposition  density  due  to  the  Implantation  is  also  shown 
(right  hand  scale).  The  arrow  marks  the  as-implanted 
amorphous-crystalline  interface  (see  text). 
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boron  profile  except  at  depths  less  than  600  X,  where  the  electrical  pro- 
file  is  larger  by  a  factor  of  2-3.  After  a  triple  sweep  at  64  W/cm  ,  the 
profile  decays  and  becomes  unmeasurable  beyond  825  X.  Also  shown  are  con¬ 
ductivity  mobilities  calculated  assuming  uc  =  0.73  pjj.  The  solid  curve 
corresponds  to  mobilities  achieved  in  bulk  material  at  the  concentration  of 
the  as  implanted  atomic  profile  [77t  Id, 102].  Only  the  triple  sweep  profile 
deviates  significantly  from  this  bulk  mobility  profile. 

Figure  4.1  includes  an  electrical  profile  measured  by  Tsai  et  al. 
[80,81]  in  low  temperature  furnace  annealed  material,  otherwise  prepared 
under  identical  conditions.  In  the  region  beyond  about  1400  X,  they  observe 
an  inactive  tail  (shaded  region) .  This  tail  was  found  to  correspond  to 
the  region  which  was  heavily  damaged,  but  not  driven  amorphous  by  the  im¬ 
plantation.  The  original  amorphous-crystalline  boundary  is  marked  by  an 
arrow  in  Figures  4. 1-4. 5.  In  the  SLEB  annealed  sample,  relatively  high 

levels  of  activation  are  measured  out  to  1900  X,  well  into  the  inactive  tail 

2 

observed  previously.  In  the  30  W/cm  sample,  we  find  boron  activation  out 

to  2600  X,  nearly  twice  as  deep  as  observed  in  550° C  furnace  annealed  material 

The  actual  distribution  of  boron  is  revealed  in  Figures  4. 2-4.4. 

Here  the  electrical  data  of  Figures  4.1  are  compared  with  the  as- implanted 

2 

and  beam  annealed  atomic  profiles  as  measured  by  SIMS.  The  30  W/cm  annealed 

material  shows  a  very  definite  skewing  while  the  electrical  profile  indicates 

poor  activation  efficiency.  Activation  becomes  nearly  complete  in  the  64 
2 

W/cm  case  (Figure  4.3),  where  the  SIMS  and  electrical  data  come  into  good 
agreement.  The  discrepancy  between  electrical  and  atomic  profiles  at  depths 
shallower  than  1100  X  is  within  combined  measurement  errors.  Note  that  this 
atomic  boron  profile  is  skewed  away  from  the  as-implanted  profile  much  less 
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Fig.  4.3.  SIMS  and  electrical  profiles  js  in  figure  4.2,  but  for 
material  processes  at  64  W/cm  .  Activation  of  boron  is 
nearly  complete  with  only  minimal  redistribution. 
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than  is  the  corresponding  profile  in  the  30  W/cm  case.  Also  note  that  the 

atomic  boron  concentration  at  depths  less  than  400  X  actually  exceeds  the 
2 

30  W/cm  profile  by  factors  of  as  much  as  2.  A  much  more  pronounced  surface 
accumulation  effect  occurs  in  the  triple  sweep  case  (Figure  4.4).  This 
accumulation  is  accompanied  by  a  large  shift  in  the  tail  of  the  distribution 
toward  the  sample  surface. 

Similar  SIMS  profiles  of  implanted  fluorine  are  presented  in 

Figure  4.5.  The  beam  annealed  fluorine  curves  possess  the  general  shapes 

of  the  corresponding  boron  SIMS  profiles,  with  one  rather  striking  exception. 
2 

The  64  W/cm  profile  shows  an  abrupt  spike  of  fluorine,  peaking  at  about 

1300  X.  Also  shown  is  a  fluorine  profile  measured  by  Tsai  et  al.  [82]  in 

material  which  was  identically  implanted  but  furnace  annealed  at  900°C  for 

30  minutes.  This  furnace  annealed  sample  has  a  fluorine  profile  with  two 

2 

peaks,  centered  about  the  peak  of  the  64  W/cm  E-beam  annealed  profile. 

4.3  Discussion 

Of  the  effects  observed  here,  the  presence  of  electrically  active 

boron  in  the  original  amorphous-crystalline  transition  region  is  perhaps 

the  most  significant.  This  can  be  understood  in  terms  of  the  study  of  f elf- 

implanted  amorphous  silicon  discussed  in  chapter  3.  The  data  show  that 

under  proper  conditions  the  defect  concentration  in  the  transition  region 

can  be  reduced  to  levels  below  that  achievable  even  by  high  temperature 

furnace  processes.  This  explains  the  ability  of  the  SLEB  process  to  activate 

2 

boron  in  this  damaged  region,  as  demonstrated  in  the  30  and  64  W/cm  cases. 
Recall  from  chapter  3  that  annealing  conditions  producing  reduced  residual 
damage  also  produced  minimal  redistribution  of  dopant  atoms.  This  accounts 
for  the  reduced  skewing  effects  observed  in  Figure  4.3  as  compared  with 
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Fig.  4.5.  Fluorine  atomic  profiles  for  as-implanted  and  SLEB  annealed 
material.  The  900°C,  30  min  furnace  annealed  profile  of 
Tsai  et  al. j  82|  is  included  for  comparison.  Note  the  large 
accumulation  of  fluorine  at  the  original  amorphous-crystal¬ 
line  transition  region  in  the  64  W/cm  case. 


Figure  4.2.  In  the  latter  case,  defects  have  remained  more  numerous,  re¬ 
sulting  in  the  depressed  electrical  activity  and  poor  mobility  observed  in 
the  region  600-1600  X.  At  even  lower  power  density  (19  W/cm^),  defects  are 

so  numerous  that  no  activity  can  be  detected  at  all. 

2 

In  the  30  W/cm  case  of  Figures  4.1  and  4.2,  it  appears  that  the 
boron  activation  efficiency  increases  significantly  in  the  region  beyond 
2200  X.  This  occurs  despite  depressed  boron  activity  in  the  regrown  region 
shallower  than  1600  X.  This  effect  can  be  explained  by  considering  the 
energy  deposition  density  profile  of  the  implantation  [65,102].  (Also 
shown  in  Figure  4.1.)  At  a  depth  of  2000  X,  the  deposited  energy  has  fallen 
by  a  factor  of  10  from  its  value  at  the  amorphous-crystalline  interface. 

One  should  therefore  expect  a  precipitous  drop  in  the  pre-annealing  defect 
density  in  this  region.  The  Increased  boron  electrical  activity  is  thus  at 
least  partially  due  to  a  local  increase  in  the  pre-annealing  crystal  quality. 

It  should  also  be  noticed  that  the  calculated  peak  of  the  deposited 

energy  profile  (and  therefore  the  peak  of  the  as-implanted  damage  profile) 

is  shallower  than  the  peak  of  the  as-implanted  atomic  boron  profile. 

(800  X  vs  1390  X) .  It  was  shown  in  chapter  3  that  in  both  furnace  and  SLEB 

annealed  silicon,  dopant  atoms  tend  to  migrate  to  areas  of  high  defect 

density.  All  of  the  SLEB  annealed  profiles  show  a  shift  toward  the  peak  of 

the  calculated  damage  profile,  the  most  prounounced  effect  occuring  in  the 
2 

30  W/cm  SIMS  profile  of  figure  4.2.  This  is  expected,  since  the  most 
residual  damage  should  occur  in  a  sample  annealed  at  lower  beam  power  density 
An  analogous  shift  in  the  furnace  annealed  profile  of  Tsai  et  al. [80,81]  is 
present  in  figure  4.1.  In  this  case,  however,  dopant  migration  is  also 
influenced  by  the  motion  of  the  amorphous-crystalline  interface  and  by 
residual  damage  in  the  original  Interface  region. 


The  study  of  self -implanted  silicon  (chapter  3)  shows  that  large 
concentrations  of  defects  and  heavy  migration  of  bulk  dopant  atoms  into 
the  as-implanted  interface  region  occur  in  multi-sweep  SLEB  annealed 
material.  This  effect  explains  the  triple  sweep  profile  (figure  4.4). 

Heavy  residual  damage  results  in  depressed  values  of  mobility  and  the 
large  near  surface  atomic  boron  concentration,  due  to  migration  of  boron 

i 

to  this  damage.  Migrating  bulk  phosphorus  atoms  can  compensate  implanted 

boron  atoms,  depressing  the  active  profile. 

Verification  of  phosphorus  migration  by  SIMS  is  difficult  because 

of  a  large  background  mass  31  signal  due  to  SiM.  The  lower  limit  of  SIMS 
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phosphorus  detection  is  estimated  to  be  3  x  10  cm  .  Therefore,  effects 
of  phosphorus  redistribution  during  annealing  may  be  studied  indirectly. 

For  example,  we  note  that  the  annealed  SIMS  boron  profile  of  figure  4.4 
falls  with  increasing  depth,  while  the  electrical  profile  climbs.  This 
alone  would  imply  a  general  decrease  in  defect  density,  which  should  be 

associated  with  an  Increase  in  mobility.  We  observe,  for  example,  that 

2  o 

the  30  W/cm  mobility  profile  improves  between  100  and  600  A,  and  beyond 

1600  ft  (figure  4.]),  where  the  electrical  and  SIMS  profiles  are  in  better 
agreement.  In  contrast,  the  triple  sweep  mobility  profile  (figure  4.1) 
falls  continuously.  If  a  significant  quantity  of  electrically  active  donors 
were  present  within  the  measured  profile,  compensation  of  the  boron 
acceptors  would  occur,  along  with  a  decrease  in  mobility,  as  observed  here. 
The  fact  that  this  effect  occurs  in  the  triple  sweep  case  is  consistent  with 
the  observation  in  chapter  3  that  donor  migration  is  most  severe  in  multi¬ 
sweep  annealed  material. 

Fluorine  redistribution  effects  in  furnace  annealed  BF+  implanted 
silicon  have  been  explained  in  terms  of  recrystallization-induced  aigration. 


gettering  by  radiation  damage,  and  thermal  outdiffusion. [82]  Since 

recrystallization  times  are  short  in  beam  annealing  processes,  one  would 

expect  only  the  latter  two  factors  to  be  significant  in  this  study.  Motion 

of  the  amorphous-crystalline  interface  influenced  the  damage  and  dopant 

profiles  of  only  the  lowest  temperature  furnace  annealed  sample  (section 

2 

3.2).  The  shifting  of  the  30  W/cm  fluorine  profile  (figure  4.5)  can  be 
explained  in  terms  of  gettering  at  the  same  residual  implantation- induced 
damage  we  believe  responsible  for  the  corresponding  boron  redistribution. 

At  64  W/cm  ,  the  lack  of  fluorine  at  depths  shallower  than  1100  ft  seems 
to  indicate  that  much  of  this  implantation-induced  damage  has  been  annealed 
away.  The  more  stable  damage  at  the  as-implanted  amorphous-crystalline 
interface  dominates  the  fluorine  migration  in  this  case.  There  is  excel¬ 
lent  agreement  between  the  location  of  the  SIMS  fluorine  peak  of  this  work 
and  that  of  the  as-implanted  interface  measured  by  Tsai  et  al. [80-82] 

(depth  1325  ft) .  These  authors  previously  accounted  for  the  two  peak  nature 
of  the  furnace  annealed  fluorine  profile  by  the  motion  of  the  interface 
region  during  recrystallization.  As  mentioned  above,  we  do  not  observe 
such  dual  peak  phenomena  in  E-beam  annealed  material.  Tsai  et  al.[65] 
have  also  demonstrated  that  fluorine  gettering  effects  are  strongest  on 

the  bulk  side  of  the  as-implanted  interface.  The  slight  skewing  of  the 
2 

64  W/cm  profile  coincides  with  their  observation. 

In  the  multi-sweep  case,  the  same  damage-aided  diffusion  effect 
responsible  for  the  large  surface  boron  concentrations  appears  to  produce 
the  observed  fluorine  distribution.  Although  fluorine  outdiffusion  appears 
likely,  the  present  studies  do  not  demonstrate  this  conclusively.  The 


63 


t 


%  * 


a 

i. ' 


$ 


c 


c 


flattedngof  the  nulti-sveep  curve  between  500  and  1CD0  £  tends  to  Indicate 
that  some  residual  as-implanted  Interface  damage  may  still  be  affecting 
the  final  fluorine  distribution. 

Finally,  it  should  be  noted  that  the  tails  of  the  fluorine 
profiles  move  progressively  closer  to  the  surface  with  increasing  E-beam 
treatment.  Furnace  annealing  studies [81]  have  shown  that  radiation  damage 
will  effectively  getter  fluorine,  even  if  it  is  as  much  as  3000  £  deeper 
than  the  peak  of  the  as-implanted  fluorine  profile.  Therefore,  it  is 
unlikely  that  any  significant  regions  of  damage  exist  deeper  than  the 
peak  of  the  as-implanted  fluorine  profile  to  the  limit  of  the  deepest 
SIMS  analysis  (7000  £) .  We  contend,  therefore,  that  even  triple  sweep 
SLEB  does  not  produce  bulk  damage  beyond  the  as-implanted  profile. 

4.4  Conclusions 

At  the  lowest  power  levels. studied  have,  electron  beam  heating 
is  sufficient  only  to  grossly  reorder  the  lattice.  At  higher  annealing 
power,  the  defect  density  is  reduced  sufficiently  to  produce  partial 
electrical  activation..  The  residual  damage  from  the  implantation  strongly 
Influences  the  final  atomic  distributions.  Under  more  intense  beam 
conditions,  this  residual  implantation  damage  is  more  completely  annealed 
and  diffusion  effects  become  less  Influential.  The  more  annealing-resistant 
damage  at  the  original  amorphous-crystalline  transition  region  strongly 
getters  Implanted  fluorine.  Finally,  under  miltiple  sweep  annealing 
conditions,  the  beam  Induces  large  defect  concentrations  in  the  region 
shallower  than  the  peak  of  the  implanted  profile.  Damage  aided  diffusion 
processes  skew  the  profiles  toward  the  surface.  The  possibility  of  some 
bulk  dopant  involvement  and  of  actual  outdif fusion  exists,  but  cannot  be 
confirmed  through  this  study. 


The  data  are  cons  latent  with  recrystallization  and  dopant  activation 
in  amorphous,  Bit  implanted  silicon  by  swept  line  electron  beam  annealing. 
Under  proper  conditions,  acceptor  profiles  with  good  electrical  activity 
and  mobility  are  observed,  and  redistribution  effects  are  minimized  by  the 
use  of  SLEB  annealing.  This  study  also  demonstrates  that  SLEB  processing 
can  produce  electrical  activity  in  the  amorphous-crystalline  transition 
region,  which  has  not  been  achieved  by  low  temperature  furnace  processes. 

The  Influence  of  residual  radiation  damage  on  the  redistribution  of  the 
implanted  species  is  consistent  with  that  observed  in  furnace  annealed 
material.  Our  results  show  that  SLEB  annealing  does  not  produce  bulk 
lattice  damage  out  to  5000  X  beyond  the  peak  of  the  as-implanted  profile. 


5.  SLEB  ANNEALING  OF  GaAs,  P 

1-x  x 

In  this  chapter  the  Investigation  of  SLEB  annealing  Is  extended 
to  Include  the  GaAs^_xPx  system.  This  technique  can  effectively  activate 
Implanted  dopants  and  the  nitrogen  lsoelectronic  trap  with  minimal  redis¬ 
tribution  of  impurities  or  implantation  related  damage.  In  the  first  half 
of  the  chapter,  annealing  effectiveness  will  be  explored  with  photo¬ 
luminescence  techniques.  The  second  half  will  be  devoted  to  characterization 
of  GaAs1_xPx  p-n  junctions. 

5.1  Study  of  SLEB  Annealed  GaAs^  Pi  N  by  Photoluminescence 

SLEB  annealing  effectiveness  was  studied  initially  by  photo¬ 
luminescence  techniques.  This  effort  was  pursued  in  cooperation  with  Dr. 

T.  Yu. [91]  The  relative  ease  of  evaluation  allowed  a  large  range  of 
annealing  parameters  to  be  explored.  The  relative  strength  of  band-edge 
and  nitrogen  trap  related  emissions  were  used  as  a  measure  of  annealing 
effectiveness.  Stripping  techniques  were  used  to  evaluate  migration  of 
nitrogen  and  implantation  related  damage. 

5.1  Experimental  procedures 

All  crystals  used  in  this  study  were  Monsanto  GaAs^_xPx  n-type 
VPE  layers  of  <100  >  orientation,  grown  on  n+  GaAs  (for  x*0.A,0.5)  or 
n+  GaP  (for  x-0.65)  substrates.  Epitaxial  layer  thicknesses  of  approxi¬ 
mately  120  pm  included  a  20-50  pm  graded  composition  region.  The  surface 
of  the  x-0.4  and  0.5  material  was  chemically  etched  to  a  depth  of  3  pm  to 
ensure  uniform  crystal  composition.  The  20  pm  just  below  the  surface  of 
the  GaASg  §5  layers  were  doPed  with  nitrogen  during  growth.  This 


doped  layer  was  etched  away  on  samples  to  be  used  In  N-  Implantation  studies. 

The  photoluminescence  emission  from  all  of  the  etched  surfaces  was  measured 

prior  to  implantation  for  reference. 

14  + 

Nitrogen  (  N  )  implantations  were  performed  at  room  temperature 

i 

at  200  keV  as  described  in  section  2.2.  Ion  doses  were  chosen  to  produce 
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peak  concentrations  of  10  -  10  cm  ,  as  determined  by  LSS  calculations. [104] 
Prior  to  electron  beam  annealing,  1000  &  of  SIO2  was  deposited  on  the  Im¬ 
planted  surface  by  the  oxidation  of  silane.  Depositions  were  typically 
performed  at  450°C  for  6  min. 

Electron  beam  annealing  was  performed  as  described  in  section  2.1. 

In  this  study,  the  beam  energy  was  maintained  at  20  keV,  and  the  table  trans¬ 
lation  speed  at  0.25  cm/ sec.  Actual  on- sample  beam  power  densities  were  not 
measured.  Rather,  total  electron  current  leaving  the  gun(I^)  was  used  as 
a  relative  measure  of  beam  strength. 

Conventional  furnace  annealing  was  performed  on  similar  samples 
for  comparison.  SiO^  films  of  thickness  1500  &  were  deposited  on  both  sides 
of  these  samples  to  prevent  decomposition  during  high-temperature  processing. 
Annealing  was  performed  in  flowing  high-purity  argon  in  a  silica-lined 
furnace.  Each  anneal  required  30  minutes  at  selected  temperatures  from 
600°  to  1000°C.  These  procedures  for  furnace  annealing  of  N- implanted 
GaAs^_xPx  have  been  discussed  in  detail  previously. [105-110] 

Prior  to  photo luminescence  measurements,  the  protective  oxides 
were  removed  with  30%  HF.  Samples  were  then  mounted  in  a  Janis  "Super 
Varitemp"  gas  exchange  liquld-*helium  cryostat  and  held  at  5  K.  The 

photoexcitation  was  provided  by  the  4880  £  line  of  a  Coherent  Radiation 

3  2 

52-G  Ar  ion  laser.  An  average  incident  power  density  of  10  W/cm  was  used. 


The  front  surface  luminescence  signal  was  focused  onto  the  slit  of  a  1/2-m 
Spex  1302  double  grating  monochromator,  and  was  detected  by  an  RCA  C31034 
(extended  S-20  response) photomulitplier  using  lock-in  amplifier  techniques. 
The  emission  spectra  are  uncorrected  for  system  response. 

Successive  layer  removal  required  for  the  profiling  study  was 
acconq>lished  by  chemical  etching  with  a  solution  of  5  H^SO^ :  1  H^O^:  1  H^O. 
The  etching  time  was  controlled  to  remove  layers  ~  600  ft  thick.  Actual 
thicknesses  removed  were  determined  by  Dektak  mechanical  stylus  measurements 
of  selective  etched  steps. 

5.1.2  Results  and  discussion 

5. 1.2.1  Direct  gap  composition,  x*0.40 
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PH'otoluninescence  spectra  for  10  cm  nitrogen  implanted  samples 
as  a  function  of  annealing  electron  beam  current  (Ib>  are  shown  in 
figure  5.1.  We  observe  as  much  as  a  30  ft  variation  of  the  near  band-edge 
emission  across  any  given  sample  owing  to  slight  variations  in  alloy  com¬ 
position.  This  accounts  for  the  small  peak  misalignments  observed  in 
spectral  data.  Shallow  donor  band-edge  emission  T  (X-  6160  ft)  becomes 
significant  after  annealing  with  Ib“10  ma  (curve  a) .  As  the  annealing 
current  is  increased,  strong  nitrogen  activation  is  indicated  by  an  increase 
in  Nr(X-6250  X) [105,106]  and  Nx(X-6500  X) [105-108]  emission.  At  Ib-15  mA, 
we  observe  localized  surface  melting  and  a  reduction  of  emission  from 
these  melted  areas.  We  observe,  however,  that  significant  increases  in 


emission  are  possible  with  two  sweeps  of  the  beam  with  Ibm12  mA.  Nx 

emission  becomes  dominant  (curve  d)  and  total  integrated  emission  becomes 

comparable  to  that  of  the  best  furnace  annealed  samples  (curve  e) .  Annealing 

19 

characteristics  of  samples  with  peak  nitrogen  concentrations  of  10  and 
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10  cm  are  similar. 


Fig  5.1.  Photoluminescence  spectra  of  GaAs.  ,Pn  .  nitrogen 

U.O  0.4  jg  __ 

implanted  to  a  peak  concentration  of  10  cm  and 
swept-line  electron  beam  annealed  at  the  indicated 
total  beam  current  (I.  ).  Emission  from  unlmplanted 
nitrogen  free  material  is  shown  for  comparison 
(dashed  line) .  (After  Yu  I  91 1 . ) 


For  peak  nitrogen  concentrations  of  10  cm,  the  maximum  emis¬ 
sion  is  observed  after  a  single  annealing  pass  at  1-12  mA.  The  N-induced 

o 

emission  intensity  increases  slightly  with  increasing  I,  or  multiple  sweep- 

b 

ing,  while  the  total  integrated  intensity  remains  the  same.  This  suggests 
that  much  of  the  lattice  disorder  introduced  by  implantation  has  been 
removed  at  1^*12  mA,  but  that  larger  1^  is  necessary  for  efficient  nitrogen 
activation. 

The  data  for  total  integrated  emission  intensity  versus  furnace 
annealing  temperature  agree  generally  with  those  observed  by  Anderson 
et  al. [109,110]  The  integrated  intensities  from  samples  electron  beam 
annealed  under  optimum  conditions  (1-12  mA,  two  passes)  are  comparable 

D 

or  superior  to  furnace  annealed  samples  at  their  individual  optimum  conditions 

This  is  true  for  all  peak  concentrations  of  implanted  nitrogen.  A  high 
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degree  of  nitrogen  substitution  is  evident  in  the  case  of  10  cm,  for  which 

total  integrated  emission  exceeds  that  of  unimplanted  direct-gap  materials. 

5. 1.2. 2  Indirect  gap  composition,  x*0.  50 

The  annealing  characteristics  of  indirect-gap  GaAs^  ^  are 

very  similar  to  those  obtained  for  direct  gap  x-0.40  material.  A  weak 

Dx  (A-5970  X,  exiton  bound  to  neutral  donor)  emission  appears  after  anneal- 
o 

ing  with  1^-10  mA.  As  the  annealing  current  is  increased,  the  emission 
(A-6350  X)  strengthens  until  it  saturates  after  a  double  anneal  at  1^-12  mA 
(curve  d) .  Integrated  intensity  obtained  under  this  annealing  condition 
is  comparable  to  that  obtained  from  optimum  furnace  annealing  (950°C,  30 


min).  This  trend  holds  for  all  nitrogen  concentrations  studied  here. 
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The  spectral  shapes  observed  are  essentially  the  same  as  furnace  annealed 
samples  for  all  nitrogen  concentrations  studied. 

5. 1.2. 3  Indirect  gap  composition,  x»0.65 

Similar  annealing  studies  were  performed  on  GaAs.  P  samples. 

U • jJ  0 • Dj 

D*  emission  (X»5810  X)  and  emission  (X-6010  X)  become  strong  under  the 

same  annealing  conditions  as  observed  In  x»0.4  and  0.5  material.  Again  we 

find  optimum  E-beam  annealing  to  produce  emission  intensity  comparable 

to  optimum  furnace  annealed  samples.  We  also  find  the  peak  Intensity  of 
19  -3 

10  cm  Implanted  material  to  exceed  that  obtained  from  similar  material 
doped  with  nitrogen  during  growth.  This  demonstrates  the  effectiveness  of 
the  electron  beam  process  in  nitrogen  activation  and  damage  annealing. 

5. 1.2. 4  Optical  depth  Profiling 

Depth  profiling  studies[91]  were  performed  on  x-0.4  and  0.5  material 
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with  peak  nitrogen  concetfcrations  of  10  cm.  The  x->0.65  material  was 
not  studied  because  of  its  unsatisfactory  etching  characteristics.  Figure 
5.2  shows  total  integrated  emission  intensity  for  both  annealed  and  unan¬ 
nealed  profiles  with  x*0.4.  Annealing  conditions  were  chosen  to  produce 
maximum  integrated  emission.  Variation  in  intensity  was  found  not  to 
exceed  20%  across  any  etched  region. 

The  unannealed  x*0.4  material  produces  no  detectable  emission 
until  the  first  20  urn  have  been  removed.  Band-edge  T  emission  is  observed 
to  increase  until  0.58  ym,  where  the  intensity  recovers  to  the  level  of  the 
unimplanted  material.  The  damage  profile  for  this  implantation  is  predict¬ 
ed  by  Brice[103]  to  have  an  asymmetrical  shape  peaking  at  2600  X.  The 
damage  remaining  at  0.58  ym  should  be  only  9%  of  the  peak  damage.  Our 
measurements  show  that  the  damage  profile  before  annealing  is  no  deeper 
than  predicted.  For  annealed  samples,  the  total  photoluminescence  emis- 


Normalized  Integrated  Intensity 
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Fig  5.2.  Depth  dependence  of  the  total  photoluminescence  intensity  (5°K) 
of  GaAs.  ,P_  ^  nitrogen  implanted  to  10  ctn“3  peak  concentration. 
Dashed  curve* corresponds  to  furnace  annealed  material  (950°  C, 

30  min),  and  solid  line  to  SLEB  annealed  material  (I.  *12  mA, 
double  pass).  Data  are  normalized  to  unimplanted,  unannealed 
nitrogen-free  material.  (After  Yu  | 91 f . > 
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sion  intensity  in  the  first  0.2  y®  is  found  to  be  comparable  with  that  of 
the  unimplanted  materials,  owing  to  the  high  concentration  of  active 
nitrogen  in  this  region.  Beyond  this  point  emission  decreases,  with  recov¬ 
ery  to  the  level  of  unimplanted  material  occuripg  at  0.6  and  1.2  ym, 
respectively,  for  the  electron  beam  and  furnace  annealed  samples.  The  same 
general  trend  is  found  in  x->0.5  material.  These  dips  are  due  to  a  combination 
of  reduced  nitrogen  concentration  and  residual  implantation-induced  damage. 

The  unannealed  sample  shows  no  significant  emission  degradation  at  depths 
of  0.6  to  1.2  ym.  The  degradation  of  the  furnace  annealed  samples  in  this 
same  region  may  therefore  be  attributed  to  diffusion  of  defects.  The  short 
duration  of  the  electron  beam  anneal  prevents  significant  damage  diffusion, 
as  demonstrated  by  the  more  rapid  recovery  of  emission  with  depth  into  the 
sample. 

In  addition  to  providing  information  about  defect  distributions, 
these  optical  profiling  measurements  can  be  used  to  study  the  effects  of 
annealing  on  the  Implanted  nitrogen  prof ile. [109]  The  photoluminescence 
intensities  obtained  after  successive  etch  steps  depend  upon  the  overlap 
of  the  optically  generated  excess  carrier  distribution  and  the  nitrogen 
profile.  The  active  impurity  profiles  may  be  deduced  approximately  if 
the  excess  carrier  profile  is  known.  Here  we  use  the  expression  of  6p(x) 
due  to  Williams  and  Chapman. [Ill]  For  x-0.4  material  we  take  o(X»4880  £)* 

4.2  x  10^  cm~^  for  5°K.  For  x«0.5,  we  take  a(4880  X)  to  be  3.3  x  10^  cm  * 
at  this  same  temperature. [91]  The  diffusion  lengths  are  expected  to  be 
very  short  compared  to  the  absorption  lengths  due  to  rapid  trapping  at  the 
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N  center  and  the  effects  of  lattice  damage.  We  use  a  radiative  lifetime 
for  the  N-trap  to  be  about  100  ns [112],  a  surface  recombination  velocity 
of  5  x  10^  cm-s  \  and  a  diffusion  length  of  ~  0.2  pm. 

The  measured  Remission  profiles  for  x*0.4  and  0.5  material  are 
shown  in  figures  5.3  and  5.4  respectively.  Also  shown  are  the  theoretical 
emission  profiles  based  on  the  overlap  of  the  estimated  excess  carrier 
distribution  and  a  simple  Gaussian  as-implanted  nitrogen  distribution 
(Rp*3700  X,  ARp*1150  X)  predicted  by  LSS  theory.  This  computed  profile  has 
a  peak  at  ~  1200  X  for  both  x»0.4  and  0.5  compositions.  There  is  good 
agreement  between  the  theorectical  active  nitrogen  profile  and  the  electron 
beam  annealed  profile.  The  emission  in  both  x*0.4  and  0.5  E-beam  annealed 
material  ceases  abruptly  after  removal  of  the  implanted  layer.  This  is 
consistent  with  the  minimal  diffusion  expected  in  solid  phase  transient  an¬ 
nealing  which  has  been  demonstrated  in  chapters  3  and  4. 

Although  active  nitrogen  profiles  for  the  furnace  annealed  samples 
reach  a  peak  at  the  same  depth  as  the  E-beam  annealed  samples,  they  display 
detectable  Nx  emission  to  depths  of  0.8  and  1.2  pm  respectively,  for 
x-0.4  and  0.5  material.  For  x«0.4  material,  0.8  um  actually  represents 
the  point  where  we  are  no  longer  able  to  discern  emission  above  the 
donor-acceptor  pair  and  Np  emission.  Since  we  observe  Np  emission  to 
a  depth  of  1.2  pm  in  the  x*0.4  case,  it  is  likely  that  nitrogen  has  dif¬ 
fused  as  deeply  as  in  the  x**0.5  composition  material.  Since  pair  emission 
is  negligible  in  the  case  of  indirect  material  and  there  is  no  Nr  inter- 
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ference,  we  can  unambiguously  define  the  limit  of  nitrogen  diffusion  with 
the  profile.  We  estimate  that  the  diffusive  tails  contain  about  10Z 
of  the  total  implanted  nitrogen  for  either  composition. 

5.2  Study  of  SLEB  Annealed  GaAsQ  55  p-n  Junctions 

The  study  described  in  the  preceding  section  shows  the  effective¬ 
ness  of  SLEB  processing  in  activating  N-implanted  GaAs  However,  £n 

the  case  of  x>0.65  composition  material,  migration  of  nitrogen  and  damage 
migration  effects  could  not  be  verified.  In  addition,  the  Nx  intensity 
profiles  discussed  above  must  rely  upon  a  calculated  excess  carrier  profile, 
rather  than  some  more  direct  means.  To  investigate  this  particular  composition 
further,  a  series  of  GaAsQ  P~n  junctions  were  fabricated  which  in¬ 

corporated  implanted  nitrogen.  The  electrical  and  optical  characteristics 
of  these  devices  show  that  SLEB  annealing  does  indeed  reduce  migration 
effects  of  implantation  related  damage  and  nitrogen  when  compared  with  con¬ 
ventional  furnace  annealing. 

5.2.1  Experimental  procedure 

Material  used  in  this  portion  of  the  study  was  indentical  to  the 
x**0.65  composition  GaAs^^P^  used  for  the  photoluminescence  experiments. 

As  before,  the  20  pm  nitrogen  doped  epitaxial  layer  was  removed  prior  to 

processing.  After  cleaning,  samples  were  subjected  to  two  beryllium  im- 

14  2  14  2 

plantations  of  1  x  10  ions/cm  at  50  keV  and  6.7  x  10  ions/ cm  at 

14  2 

130  keV.  A  single  3  x  10  ions/cm  ,  220  keV  nitrogen  implantation  was 
then  made.  The  calculated  as-implanted  distributions  are  shown  in  figure 
5.5.  The  Be  doses  and  energies  were  chosen  to  produce  the  retrograde 
acceptor  profile  studied  by  Chatterjee. [113]  This  distribution  was  a 
compromise  between  a  uniform  and  buried  gausslan  profile.  The  former 


.a  .  * 


I 


As  -  Implanted  Profiles 
GaASQ  35PQ  53:  Be  +  N  Diodes 

Be:  50  keV,  lxl014cm’2 
130  keV,  6.7xl014cm“2 
N:  220keV,  3x  I014cm”2 


\\ 


_ Nitrogen 

_ Beryllium 


.Sulfur  (Background) 


DEPTH  (/xm) 


Fig.  5.5.  Calculated  as-implanted  profiles  of  nitrogen  and  beryllium 
for  the  GaAsQ  ^5^0  gg  p-n  junctions.  Also  shown  is  the 

level  of  the  substrate  dopant  concentration. 


distribution  simulates  the  zinc  profile  created  in  conventional  furnace 

diffusion,  whereas  the  latter  produces  a  built-in  retarding  field  for  injected 

electrons,  reducing  surface  recombination  effects.  Chatterjee’s  studies 

have  suggested  that  the  retrograde  acceptor  profile  does  indeed  reduce  surface 

recombination  and  results  in  the  best  light  intensity  (L)  versus  current 

characteristics  among  similar  implanted  direct-gap  GaAsP  light  emitting 
« 

diodes.  The  as-implanted  nitrogen  profile  in  figure  5. 5  is  entirely 
contained  within  the  implanted  acceptor  profile. 

Following  implantation,  samples  were  cleaned  and  encapsulated 
with  ~  1500  X  of  rf  plasma  deposited  silicon  nitride. £114]  Some  of  the 
material  was  then  annealed  in  silica  lined  furnace  tube  under  flowing  high 
purity  argon  gas.  A  30  minute,  950°C  annealing  cycle  was  selected,  cor¬ 
responding  to  the  optimal  conditions  for  annealing  nitrogen  implanted  to  this 
level. [91]  Chat terjee[ 113]  demonstrated  successful  annealing  of  Be- 
implanted  direct  gap  GaAsP  with  these  approximate  implantation  conditions 
at  temperatures  about  50°C  less  than  used  here.  Samples  for  E-beam  anneal¬ 
ing  were  prepared  as  described  in  section  2.1.  A  table  translation  speed 
of  0.60  cm/ sec  was  used  for  these  samples.  Although  a  number  of  different 

beam  power  density  combinations  were  tried,  the  devices  described  here  were 

2 

processed  with  five  passes  of  the  SLEB  with  18  W/cm  neak  power  density. 

Prior  to  annealing,  all  implanted  samples  possessed  a  frosted 
opaque  appearance.  After  the  first  pass  of  the  SLEB,  the  samples  returned 
to  their  original  orange  transparent  condition.  This  observation  tends  to 
indicate  that  good  crystallinity  had  returned  after  this  level  of  annealing. 

A  survey  of  photoluminescence  of  this  Be  +  N  implanted  material  agreed 
with  findings  of  section  5.1  In  general,  restoration  of  crystal  order 


occured  after  only  a  single  pass  of  the  SLEB,  but  complete  nitrogen 
activation  required  more  intense  and  repetitive  E-beam  . treatment. 


After  annealing,  samples  were  solvent  cleaned  and  stripped  of 
nitrides.  A1  contacts  were  then  defined  on  the  implanted  surface  by  AZ 
photoresist  lift-off  technique.  On  the  back  surface,  Au-Ge  contacts 
were  produced  with  Ni  overlay  by  double  evaporation  through  a  shadow 
mask.  The  sauries  were  then  sintered  in  flowing  for  a  total  of  ap¬ 
proximately  90  secs.  Maximum  temperature  reached  during  the  contact 
processing  did  not  exceed  400°C.  Individual  devices  were  formed  by  scribing 
and  cleaving.  Devices  were  mounted  in  pairs  on  T0-18  headers  and 
thermocompression  lead  bonded. 

5.2.2  Results  and  discussion;  etlectrical  characteristics 

Current-voltage  measurements  were  made  on  these  devices  with  a 

computer  controlled  Hewlett  Packard  69321  B  power  supply  and  a  Keithley 

Model  480  Picoammeter.  Temperature  dependent  measurements  were  made  in  a 

vacuum  chamber  equipped  with  an  Air  Products  closed-cycle  He  refrigerator 

unit.  Room  temperature  current  density  (J)  versus  voltage  characteristics 

of  the  best  furnace  and  SLEB  annealed  devices  are  shown  in  figure  5.6. 

In  the  low  forward  bias  region,  both  devices  display  conduction 

with  an  ideality  factor  n  of  2.*  The  furnace  annealed  device  shows  a 

-5  2 

marked  departure  from  this  behavior  at  J*10  A/ cm  ,  and  yet  further  bend- 
over  occuring  at  forward  bias  exceeding  2.5  volts.  This  departure  from 


£ 

n  is  the  ideality  factor  commonly  used  in  the  diode  expression 

J  -  J  exp(qV/nkt) 
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Fig.  5.6.  Current  density  versus  voltage  characteristics  for  the  best 
GaASg  jjPq  65:®e  +  ^  diodes  annealed  by  SLEB  and  furnace 

techniques. 
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n*2  conducCion  occurs  consistently  for  all  5  furnace  processed  devices 

-5  -6  2 

between  10  and  10  A/cm  .  The  e-beam  device  characteristic  shows 
n*2  conduction  nearly  to  the  limit  of  current  measurement.  Some  slight 
roll-off  is  seen  at  the  highest  measured  currents,  presumably  due  to 
series  resistance  effects.  All  but  one  of  the  eight  beam  processed 
devices  show  a  predominance  of  current  due  to  recombination  in  the  space 
charge  layer.  The  characteristics  of  the  furnace  annealed  devices 
Indicate  the  presence  of  series  resistance  in  the  junction. 

Eliminating  the  worst  device  from  each  set,  the  average  V^-10  V 

-7  2 

reverse  current  density  is  about  4  x  10  A/ cm  in  the  furnace  annealed 

—8  2 

devices  and  6  x  10  A/ cm  in  the  SLEB  case.  In  general,  the  furnace 

annealed  devices  show  a  softer  reverse  bias  characteristic. 

Temperature  dependent  J-V  data  for  devices  T4B  (f urance)  and  ElB 
(SLEB)  are  presented  in  figures  5.7  and  5.8  respectively.  Notice  in 
figure  5. 7  a  set  of  parallel  lines  has  been  drawn  through  the  approximate 
straight  lint-  portion  of  each  curve,  beginning  just  above  the  "knee" 
discussed  earlier.  The  obvious  temperature  independence  indicates  conduction 
predominated  by  tunneling  mechanisms. [115,116]  One  would  therefore  expect 
a  relatively  large  number  of  states  to  be  available  within  the  bandgap  at 
the  n-type  edge  of  the  space  charge  region  in  these  devices.  The  doping 
level  of  the  substrate  material  is  not  sufficiently  large  to  expect  valence 
bandtail  effects. [117]  A  similar  set  of  parallel  lines  obviously  cannot  be 
drawn  through  the  data  of  device  ElB.  Calculated  ideality  factors  of  the 
beam  annealed  device  vary  around  2  for  all  temperatures  studied. 
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Fig.  5.7.  Temperature  dependent  current  density  versus  voltage 

characteristics  for  furnace  annealed  GaAs^  ggSBe+N 

diode  #T4B.  The  solid  lines  drawn  through  the  low  current 
density  portion  of  each  curve  are  exactly  parallel  to  one 
another. 
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Fig.  S.8.  Temperature  dependent  current  density  versus  voltage 

characteristics  for  SLGB  annealed  device  GIB.  The  solid 
lines  dravn  through  each  characteristic  have  ideality 
factors  of  about  2  at  all  temperatures.  Note  the  shift  in 
horizontal  scale  from  figure  5.7. 


DLTS  analysis  tends  to  confirm  these  findings.  The  two-diode 

method  described  in  section  2.3.2  is  used.  A  survey  of  the  electron  trap 

spectra  on  the  n-type  side  of  the  junctions  of  two  representative  devices 

is  shown  in  figure  ,5.9.  Here  the  SLEB  and  furnace  annealed  DLTS  signal  are 

plotted  logarithmically  in  the  same  units.  Quiescent  reverse  bias  of  -10  V 

and  bias  reducing  pulses  of  +10  V  were  used  for  both  devices.  Corrections 

for  device  area  are  included  so  that  the  relative  magnitudes  of  these  curves 

are  indicative  of  the  relative  trap  densities  between  devices.  Also  plotted 

are  survey  spectra  of  Day [71]  measured  un  A1  Schottky  barriers  fabricated  on 

material  identical  to  that  used  in  this  study.  These  data  were  measured 

under  identical  frequency  conditions,  but  not  necessarily  the  same  bias. 

Therefore  for  comparison,  Day's  data  have  been  normalized  to  the  peak  value 

of  the  beam  annealed  sample's  spectra. 

The  SLEB  annealed  sample  shows  a  single  large  peak,  while  the 

furnace  annealed  device  shows  significant  signal  all  across  the  scan. 

The  defect  spectra  of  the  unlmplated  "nitrogen  free"  Schottky  diode  shows 

an  amazing  similarity  to  that  of  the  beam  annealed  device.  This  defect 

(Ea,,  E-E  -0.41  eV)  has  been  associated  with  the  dopant  sulfur. [118]  I 
oi  c 

conclude  that  these  defects  are  intrinsic  to  the  material  and  are  not  a 
product  of  the  annealing.  The  average  concentration  of  this  trap  can  be 
estimated  through  [69] 

N  -2  • (N  ,  -  N  ) •  AC/ C  (5.1) 

i  a  a 
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to  be  about  1  x  10  cm 


Fig  5.9.  DLTS  spectra  for  GaAs.  oc^q  gs5®6***  diodes  processed  by  SLEB  and  furnace 

annealing.  Also  shownare  normalized  spectra  for  unlmplanted  and  N-lmplanted 
x-0.65  GaAs,  P  Schottky  barrier  diodes  as  measured  by  Day  et  al.  | 7l| .  Noti 
the  logarit male  vertical  scale. 
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The  influence  of  crap  signal  Eg^  is  apparent  in  the  furnace 

annealed  *e  +  N  diode  spectra  at  about  -95°C.  This  device  also  shows  the 

effects  of  levels  and  observed  in  a  nitrogen  Implanted  Schottky 

diode. [71]  Using  equation  5.1  we  may  estimate  the  concentration  of  these 
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defects  to  be  about  2  x  10  cm  .  Although  identification  is  not 

absolutely  conclusive,  it  is  probable  that  these  signals  are  due  to  implanted 

% 

nitrogen  which  as  diffused  out  of  the  original  Implanted  profile  and  beyond 

the  metallurgical  junction  during  furnace  annealing. 

Under  similar  implantation  and  furnace  annealing  conditions, 

Chatterjee[113]  found  evidence  of  Be  diffusion  beyond  the  as- implanted 

distribution.  This  has  probably  occured  in  the  furnace  annealed  devices 

under  study  here.  They  show  consistently  high  voltage  drops  at  large 

-1  2 

forward  currents.  The  average  applied  bias  required  for  10  A/ cm  current 
density  is  3.1  V  for  furnace  annealed  devices  and  only  1.5  V  for  SLEB 
annealed  devices.  This  tends  to  indicate  a  higher  sheet  resistance  due 
to  spreading  of  the  Implanted  profile.  The  nitrogen  diffusion  under  con¬ 
sideration  here  must  have  proceeded  at  a  faster  rate  than  that  of  beryllium, 
since  it  is  detected  in  the  space  charge  beyond  the  metallurgical  junction. 

It  is  significant  that  no  nitrogen  related  signal  is  detected  in  the  SLEB 
annealed  device.  This  coincides  with  the  finding  of  the  photoluminescence 
study  of  section  5.1. 

Also  observed  is  a  broad  signal  band  in  the  furnace  annealed 

Be  +  N  diode  spectra  at  temperatures  above  -50°C.  Day [71]  has  found 

similar  broad  defect  peaks  in  nitrogen  implanted  x-0.65  GaAs^_xPx.  This 
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is  shown  in  the  truncated  5  x  10  cm  nitrogen  implanted  spectra  cf 

18 

figure  5.9.  (This  truncated  curve  is  identical  to  the  10  N-implanted 
spectra  below  -80°C.)  It  is  significant  to  note  that  the  SLEB  annealed 
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device  shows  no  similar  radiation  damage  signal  spectra,  even  though  the 
material  was  implanted  to  a  factor  of  two  larger  peak  nitrogen  concentration 
and  was  beryllium  implanted  as  well.  However,  we  cannot  be  completely 
certain  that  there  are  no  hole  traps  in ''the  depletion  layer  since  hole 
injection  was  not  performed  during  DLTS  analysis. 

Along  this  same  line  of  argument,  we  may  not  conclude  that  the 

defects  observed  here  are  directly  responsible  for  the  modes  of  current 

conduction  discussed  above.  Holes  injected  into  the  n-side  of  these 

devices  will  dominate  conduction.  However,  it  is  Interesting  that,  these 

furnace  annealed  Be  +  N  diodes,  which  possess  a  continuum  of  deep  levels, 

also  demonstrate  tunneling  current  components.  As  an  order  of  magnitude 

estimate,  if  a  200  X  thick  layer  contained  an  average  of  1  x  1015  traps/cm^ 

each  would  need  to  recombine  a  hole-electron  pai~  only  once  every  300 
_6  2 

yses  to  support  10  A/ cm  conduction.  This  is  certainly  not  an  unreasonably 

short  lifetime. 

5.2.3  Results  and  discussion:  light  emission  characteristics 


Although  the  principal  purpose  of  these  diodes  was  the  study  of 
redistribution  phenomena  by  electrical  techniques,  their  light  emission 
characteristics  also  lend  some  support  to  the  contentions  of  the  preceding 
section.  One  would  not  expect  large  emission  efficiencies  in  these  devices. 
In  this  composition,  GaAsP  has  an  indirect  gap,  making  bandedge  radiative 
emission  a  low  probability  event.  Nitrogen  acts  as  an  efficient 
isoelectronic  trap  of  this  composition,  but  the  doping  profile  has  not 


been  optimized  for  maximum  efficiency  in  these  devices. 
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Light  intensity  (L)  versus  current  data  were  taken  by  positioning 
the  devices  in  a  180°  polished  aluminum  reflector.  This  reflector  was  then 
attached  to  a  Newport  Model  880  radiometric  detector  system.  The  reflector/ 
detector  combination  was  designed  to  eliminate  all  background  optical 
signals.  Total  optical  power  measurements  were  then  made  as  a  function  of 
forward  DC  current.  The  data  are  corrected  for  photodetector  response. 

Two  representative  device  characteristics  are  shown  in  figure  5.10.  Shown 
for  reference  are  lines  indicating  emission  proportional  to  current  density 
and  to  current  density  squared. 

The  slopes  of  the  characteristics  shown  are  representative  of  all 
devices  annealed  by  the  respective  techniques.  In  general,  furnace  annealed 
diodes  all  show  behavior  while  all  SLEB  annealed  devices  show 

jf\»l 

L°cJ  .  This  trend  Indicates  that  some  space  charge  recombination  current 
is  contributing  to  the  emission  of  the  furnace  annealed  diodes,  while  dif¬ 
fusion  current  related  emission  dominates  in  the  SLEB  annealed  case. [113] 

This  finding  is  not  Inconsistent  with  the  J-V  measurements  discussed 
earlier.  Ohmic  losses  are  considerable  at  the  current  densities  where 
light  intensity  is  large  enough  to  measure.  Therefore,  we  are  not  able  to 
identify  the  dominant  mode  of  conduction  at  these  current  levels.  The 
average  external  quantum  efficiency  for  the  furnace  annealed  devices  is 
about  5  x  10  ®  photons/electron  for  the  furnace  annealed  diodes  and  about 
2.8  x  10”^  for  SLEB  annealed  devices.  A  figure  of  about  2  x  10  5  is  expected 
from  published  results  in  furnace  diffused  x*0. 65  LEDs  without  nitrogen 
doping. [119] 
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Since  nitrogen  is  known  to  exist  in  the  space  charge  region  of  the 
furnace  treated  devices,  one  might  expect  that  some  isoelectronic  trap 
emission  might  be  observed.  This  occurs  to  some  extent.  The  detection 
system  described  in  section  5.1  was  used  to  make  electroluminescence 
measurements.  The  peak  emission  wavelength  of  the  beam  annealed  device 
E1A  6018  X  (2.060  eV) ,  is  consistent  with  that  expected  from  measurements 

I 

of  Crawford  et  al. [120]  for  GaAsP  without  nitrogen  doping.  Very  little 

current  dependent  emission  shift  is  observed  in  this  device  from  17  to 
2 

40  A/cm  .  Peak  emission  for  the  furnace  annealed  device  T4A  is  nearly 

identical  to  E1A  at  low  bias,  but  shifts  over  70  £  to  6078  X  (2.040  eV) 

2 

at  25  A/ cm  .  This  wavelength  corresponds  roughly  with  the  expected 
(A-line)  nitrogen  emission. [120] 

The  results  sugggest  again  that  some  depletion  region  nitrogen 
Involvement  occurs  in  the  case  of  the  furnace  annealed  devices.  It  is 
possible  that  the  nitrogen  related  emission  also  has  some  component  from 
electrons  injected  into  the  p-side  of  the  junction.  The  reduced  doping 
gradient  expected  in  the  furnace  annealed  case  would  allow  increased 
electron  injection  efficiency.  The  current  dependence  of  the  emission 
tends  to  confirm  this  proposition. 

5.3  Summary  and  Conclusions 

In  this  chapter  it  was  demonstrated  that  SLEB  annealing  is  ef¬ 
fective  in  activating  ion-implanted  nitrogen  in  both  direct  and  indirect 
gap  composition  GaAs^^P^.  Photoluminescence  analysis  shows  that  under 
some  conditions  nitrogen  related  emission  can  actually  exceed  that 
attainable  by  furnace  annealing.  Photoluminescence  profiling  demonstrates 
that  SLEB  annealing  isalso  effectivein  restricting  the  migration  of  nitrogen 
and  ion-implantation  related  damage  to  the  original  implanted  layer. 
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Although  this  technique  could  not  be  used  to  verify  this  finding  in  x-0.65 
composition  material,  electrical  studies  of  Be  +  N  implanted  diodes  do  tend 
to  confirm  this  finding.  J-V,  DLTS,  L-J,  and  spectral  emission  characteristics 
all  can  be  understood  in  terms  of  restricted  motion  of  nitrogen  and  related 
damage  during  annealing.  These  data  also  show  a  lack  of  Be  migration  in 
SLEB  treated  devices,  in  contrast  with  the  furnace  annealed  case. 


CHAPTER  6  -  SUMMARY  AND  RECOMMENDATIONS 

In  this  study,  the  behavior  of  swept-line  electron  beam  (SLEB) 
annealed,  ion-implanted  semiconductors  has  been  studied  in  some  detail. 

In  general  it  is  found  that  SLEB  processing  can  be  an  effective  and 
practical  alternative  to  furnace  annealing.  SLEB  treatment  has  been 
shown  to  be  effective  in  annealing  implantation  related  damage  in 
silicon,  especially  in  the  as-implanted  amorphous-crystalline  transition 
region.  This  transition  region  damage  has  been  shown  to  be  resistant 
even  to  high  temperature  furnace  annealing,  and  can  act  to  skew  dopant 
distributions  through  damage  aided  diffusion.  These  migration  effects 
are  much  less  pronounced  in  SLEB  annealed  amorphous  silicon.  SLEB 
induced  activation  of  boron  in  this  transition  region  has  also  been 
demonstrated.  This  has  not  been  observed  in  furnace  annealed  silicon. 

The  success  of  this  annealing  technique  is  not  limited  to  silicon. 
Photoluminescence  studies  have  shown  the  effectiveness  of  the  SLEB 
technique  in  activating  Implanted  nitrogen  in  GaAs^^P^.  In  some  cases, 
nitrogen  isoelectronic  trap  related  emission  intensities  are  actually 
larger  than  those  observed  from  similarly  prepared  furnace  annealed 
material.  Photoluminescence  profiling  and  p-n  junction  studies  show  that 
SLEB  annealed  material  does  not  display  dopant  and  implantation  damage 
migration  effects  present  in  optimally  prepared  furnace  annealed  material. 

Although  this  study  demonstrates  that  SLEB  annealing  is  an  ef¬ 
fective  processing  technique,  a  number  of  questions  remain  to  be 
answered.  The  possibility  of  upscaling  a  line  beam  to  treat  industrial- 
size  wafers  was  not  completely  investigated.  A  prototype  electron  gun 
with  a  four  inch  beam  was  built  but  could  not  be  completely  developed. 

As  designed,  however,  it  meets  the  criteria  observed  to  be  necessary 


for  successful  large  scale  applications.  First,  it  does  not  rely  upon 
overlap  of  annealed  areas,  but  rather  treats  large  substrates  in  a  single 
pass.  Experiments  with  the  electron  gun  described  in  chapter  2  show  that 
stress  marks  develop  on  material  adjacent  to  the  beam  track  on  wafers 
larger  than  the  beam  spot.  No  such  stressing  appear  within  the  beam 
heated  regions.  Also,  the  long-line  gun  will  deliver  all  of  its  energy 
in  a  thin  (10  mil)  line.  This  should  improve  annealing  uniformity, 
especially  on  samples  which  are  non-rectangular.  If  the  sample  is  heated 
only  within  a  thin  line  and  is  translated  fast  enough,  significant  lateral 
heat  flow  will  not  occur.  This  should  limit  hot  spots  on  the  substrates 
which  occur  when  broader  beam  shapes  are  used. 

The  results  discussed  here  represent  only  a  relatively  limited 
portion  of  the  range  of  annealing  parameters  possible.  The  fact  that  an¬ 
nealing  quality  can  exceed  that  achievable  by  conventional  annealing 
should  stimulate  further  study.  The  question  of  the  exact  process  by 
which  point  defects  are  annealed  has  not  been  adequately  answered.  The 
degree  of  permanence  of  this  annealing  is  unknown.  Electric  field 
dependent  aging  studies  of  SLEB  annealed  devices  would  be  a  useful 
investigation.  Some  obvious  differences  in  the  annealing  conditions 
optimal  for  different  semiconductors  were  noted  during  these  studies. 

For  example,  multi-pass  annealing  depresses  the  active  boron  profiles  in 
Implanted  Si,  while  it  works  quite  well  in  implanted  GaAsP.  The  electron 
range  and  implantation  end  damage  distributions  all  undoubtedly  affect 
the  ultimate  quality  of  treatment.  Some  theoretical  treatment  of  beam- 
material  interaction  could  lead  to  more  predictable  processing. 
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TRAPSI  PROGRAM  DOCUMENTATION 

TRAPS I  is  a  FORTRAN  program  for  computing  defect  depth  profiles 
from  DLTS  measurements  at  fixed  reverse  bias  and  varying  bias  reducing 
pulse  height.  This  program  is  designed  for  analysis  of  Schottky  barriers 
fabricated  on  n-type  silicon.  Modifications  can  easily  be  made  for  other 
semiconductor  materials.  The  general  operating  principles  of  TRAPSI  were 
described  in  section  2. 3. 2. 2.  The  appendix  includes  specific  input 
requirements,  a  block  by  block  description  of  operation,  a  copy  of  the 
program,  and  a  sample  output. 

A.l  Input  Requirements 

The  TRAPSI  program  requires  three  data  types:  (i)  that  which 
relates  to  die  device  in  general;  (11)  C(V)  data  for  determination  of 
fixed  charge  concentration;  and  (iii)  the  DLTS  data  relative  to  the  defect 
level  in  question.  Data  may  be  entered  indefinitely  so  that  defect  peaks 
of  different  activation  energy  and  from  different  devices  may  be  analyzed 
without  restarting  the  program. 

Immediately  after  execution,  the  program  asks  for  type  i  data 
through  the  interactive  terminal  (#5).  The  data  are:  (The  FORTRAN 
variable  is  listed  after  each  entry.) 

a.  Descriptive  one  line  comment. 

2 

b.  Device  area  in  cm  (AREA). 

c.  Relative  dielectric  constant  for  the  semiconductor  (ER). 

The  static  value  may  be  used  with  confidence  for  Si,  Ge  and  GaAs.  Other 
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values  may  be  more  applicable  for  materials  with  relatively  large 
dielectric  relaxation  times  [121] 

d.  Metal-semiconductor  barrier  height  in  volts  (DPHIBN) . 

e.  Density  of  states  effective  mass  for  electrons  (MDEE)  and 
holes  (MDEH) . 

f.  Donor  activation  energy  in  electron  volts.  (EDON) . 

The  program  then  asks  for  the  file  name  which  contains  the  C(V) 
data  (type  ii)  for  the  device  in  question.  This  file  must  contain  two 
columns  of  data,  the  left  hand  containing  the  capacitance  in  pf,  the  right 
hand  the  corresponding  voltage  in  volts.  Up  to  100  data  pairs  may  be 
entered.  The  last  pair  of  data  points  in  the  file  must  both  be  zero. 

All  negative  voltages  should  be  entered  as  negative  numbers.  All  ef¬ 
fective  bias  voltages  and  corresponding  steady  state  capacitances  used 
in  the  DLTS  profiles  should  be  included  in  this  list.  All  static  capaci¬ 
tance  data  are  automatically  reduced  by  0.68  pf,  the  average  parallel 
capacitance  of  a  TO-18  header. 

Since  the  fixed  charge  concentration  is  inversely  proportional  to 
dC/dV,  it  is  important  that  enough  static  capacitance  data  be  entered  to 
produce  good  estimates  of  N^Cx) .  The  program  has  an  internal  criterion 
that  dC/dV  be  calculated  from  capacitance  values  which  differ  by  at  least 
IX  and  no  more  than  10%  of  the  capacitance  in  question.  If  this  criterion 
is  not  met,  an  error  flag  appears  on  the  far  left  side  of  the  output.  C(V) 
data  may  be  entered  either  largest  or  smallest  value  first,  but  must 


increase  or  decrease  continuously. 

The  defect  peak  related  data  (ty'e  iii)  are  then  entered.  They 


a.  Descriptive  comment  relative  to  this  defect. 

b.  Value  of  fixed  DC  reverse  bias  maintained  during  profiling, 
in  volts.  (VBIAS) . 

c.  The  temperature  in  °K  at  which  the  peak  trap  signal  occurs 

(TEMP). 

d.  Total  number  of  DLTS  data  points  (CPTS). 

e.  The  value  of  the  bias  reducing  puls  height  in  volts  (VR) , 
the  peak  recorder  displacement  in  cm  (DISP(I)),  the  sensitivity  of  the 
lock-in  amplifier  in  volts  (S(I)),  and  the  calibration  factor  in  V/pf 
(CVF(I)). 

The  program  will  ask  for  data  of  type  e  repetitively  until  the 
total  number  of  transient  points  (CPTS)  have  been  recieved.  After  the 
last  entry,  the  program  will  ask  if  another  defect  peak  is  to  be  analyzed. 
If  so,  it  will  ask  for  new  type  iii  data.  If  not,  another  device  may  be 
analyzed  by  supplying  new  type  i  and  ii  data.  A  sample  interactive  out¬ 
put  is  shown  on  the  following  two  pages.  Entries  made  by  the  operator 
are  underlined. 

A.  2  Block  Description 

A  description  of  TSAPSI  code  is  presented  in  this  section.  The 
line  numbers  found  at  the  beginning  of  each  entry  correspond  to  these 
generated  by  the  DEC  FORTRAN  editor  as  shown  on  the  program  listing  found 
in  section  A. 3. 

Line  100-900:  Descriptive  comments. 

Line  1000-1700:  Variable,  array  and  type  declaration. 

Line  1800-5200:  Entry  of  type  i  data. 

Line  5300-6900:  Entry  of  file  name  containing  C(V)  data  arrays 
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EXECUTE  TRAPSI.FOfi 
LINK!  Lo«tiinS 
CLNKXCT  TRAPSI  execution:) 
ENTER  ONE  LINE  COMMENT 


PROFILE  FOR  E  BEAN  S,  8.  «31A 
ENTER  DEVICE  AREA  (CM**2)  ~ 1.32E-3 

ENTER  RELATIVE  DIELECTRIC  CONST-!  ~11.8 
ENTER  DELTA  PHI  B-N  (V0LTS)=0.7 


DENSITY  OF  STATES  EFFECTIVE  fitSsES  OF  COHHON 
SEMICONDUCTORS! 

ELECTRONS 
1  .1 
.068 
.089 
1.20 
1.20 


SI 

GAA3 

GAASP  40% 

GAASP  INDIRECT 
GAP 


HOLES 

.55 

.5 


MDE  FOR  ELECTRONS -1 . 1 
MDE  FOR  HOLES- .55 

DONOR  ACTIVATION  ENERGY  <£V)-0.044 
ENTER  NAME  OF  C(V)  DATA  FILE! TFETlA 
C(V)  DATA  READ  A-OK 


ENTER  1  LINE  COHMENT 


PEAK  I  —  0.32  EV 

VALUE  OF  DC  BIAS  < VOLTS )-- 10 . 00 

ENTER  TEMP  OF  PEAK < DEG  K>  183.50 

ENTER  DELTA  CAP  DATA.  *  OFTTS^T 

BIAS  REDUCTING  PULSE(VOLTS)  V(  1)*10.5 

ENTER  VAL  *  1  OF  DISP  (CM).  SENS.  (V), 


CONVER.  FACT.  (V/PF) 


8.02  .020  t5A7 

VALUE  #  10F  PULSE=  10.30000 

OF  DISP*  8.02000 

OF  SENS*  0.02000 

OF  CON.  FACT.*  0.36700 

OK?  —  1=YES.0»N0:  1 
BIAS  REDUCTING  PULSE(VOLTS)  V<  2 ) - 1 0 ■ 00 
FNTER  VAL  *  2  OF  DISP  (CM),  SENS.  (V).  CONVER.  FACT. 


(V/PF) 


8.19  .020  .530 

VALUE  *  20F  PULSE*  10.00000 

OF  DISP=  8.19000 

OF  SENS*  0.02000 

OF  CON.  FACT.*  0.33000 

OK?  --  1*VES»0=N0:  i 

7  0.500  1 . 23290E+02  6.99840E-02 

2 . 020 1 3E-06  2 . 29372E+ 1 7  3.71246EF12 

BIAS  REDUCTING  PULSE(VOLTS)  V(  3)-A.50 
ENTER  VAL  ♦  3  OF  DISP  (CM),  SENS.  (V),  CONVER.  FACT. 


94309E-02 


(V/PF) 


6.50000 


8.27  .020  .542 

VALUE  *  ~ 37TF  PUtTST? 

OF  P I S  P  =  3.27000 

OF  SEN 8*  0.02000 

OF  CON.  FACT.*  0.34200 

OK?  --  t*YES,0*HO:  ± 

F:  •  AS  REDUCTING  FULSE(VOLTS)  V(  3>*9  ■  5 
ENTER  VAI.  *  3  OF  DISP  (CM).  SENS.  (V), 


CONVER.  FACT.  (V/PF) 
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%1 
..  * 


8.27  .020 

VALUE  *  30F  PULSE3  9.50000 

OF  DISP-  8.27000 

OF  SENS-  0.02000 

OF  CON.  FACT , =  0.00000 

OK?  —  1=YES>0=N0: 

B  t  AS  REDACTING  FULSE(VOLTS)  V<  3  )  3  9 . 5 

ENTER  VAL  *  3  OF  DISP  <CM>»  SENS.  (V).  CQNVER.  FACT.  <V/PF> 


te 


8.27  .020  . 542 

VALUE  *  30F  PULSE-  9.00000 

OF  DISP=  8.27000 

OF  SENS-  0.02000 

OF  CON.  FACT . =  0.54200 

OK?  —  1  =  YES.O-NO:  1 

2  0.000  1.00880E+02  9.24770E-02  3.84820E-02  -2.63722E-04 

9 • 67390F-06  4.28362E+16  6.76468EE12 

BIAS  REDUCTING  PULSE(VOLTS)  V<  4)39.0 

ENTER  VAL  *  4  OF  DISP  < CM ) .  SENS.  <V>#  CONVER.  FACT.  <V/'PF> 


| 


m 

-  -v' 


8.19  .070 

VALUE  •  40F  PULSE-  9.00000 

OF  DISP=  8.19000 

OF  SENS3  0.02000 

OF  CON.  FACT . -  0.36300 

0K9  —  1  -  YES  t  03N0  I  _1 

0  -0.300  7.32700E+01  9.36192E-02  4.34221E-07  2.28399E-04 

1.51 893E-03  3 • 92206E+ 16  1.31136E+13 

0  -1.000  6 . 04500E+01  9.36192E-02  4.69120E-02  -3.43861E-03 

1 . 94826E-05  3.92378E+16  3.27133E+14 

IS  THERE  ANOTHER  TRAP  PEAK  FOR  THIS  DEVICE? 
t  =  YES  *  O-NO 


_P 

Till  YOU  UISH  ro  ANALYSE  ANOTHER  DEVICE? 
1=YESf  0-NO 


END  OF  EXECUTION 

:pu  tine;  i.?6  elafsed  time:  4:21.17 

EXIT 


1 s 


(CCV(I)  and  VCV(I)).  Included  are  checks  for  end  of  file  (line  6000), 
correction  for  parallel  header  capacitance  (6300)  and  file  overflow 
(6600) . 

Line  7000-7900:  Entry  of  DLTS  defect  peak  related  comment. 

Line  8000-9100:  Entry  of  quiescent  DC  reverse  bias  (VB1AS)  and 
search  of  corresponding  value  of  static  capacitance  (CBIAS).  If  the 
search  fails  to  find  a  VCV(I)*VBIAS,  then  CCV(2)-CBIAS. 

Line  9200-10000:  Entry  of  defect  peak  temperature  (TEMP)  and 
number  of  transient  capacitance  points.  N^,  Nv>  and  E^  are  calculated 
in  the  block. 

Line  10100-10400:  Initialize  profile  calculation  accumulators. 

Line  10500:  Initialize  the  main  DO  loop  for  profile  calculation. 

i  Line  10600-10700:  Reset  dC/dV  flags.  (F,  FI) 

Line  19800:  Test  for  I*CPTS+1.  If  true,  programs  skips  over 
transient  capacitance  data  input  steps.  Note  that  DLTS  data  are  received 

I  only  when  1<I<CPTS,  and  values  of  trap  concentration  are  calculated  only 

when  2<_  I<CPTS+1. 

Line  10900-12300:  Input  and  confirm  values  for  the  type  iiie. 

!  (VR,  DISP(I) ,  S(I) ,  CVF(I) ) .  Largest  values  of  bias  reducing  pulse  (VR) 

must  be  entered  first.  Allowance  is  made  to  correct  data  if  miskeyed. 
Effective  bias  (VDELC(I))  is  calculated  in  line  11200. 

Line  12400:  Calculate  Ith  value  of  AC/C  (DELCT(I)) . 

Line  12500:  If  1-1,  VDLEC(0)-VBIAS .  This  step  allows  calculation 
of  the  fixed  charge  concentration  (NB(1))  at  the  depletion  depth  cor- 

I 

'  responding  to  VBIAS. 


Line  21600-13600:  Search  for  the  value  of  VCV(J)  which  equals 
VDELC(I-l) .  If  the  search  fails,  the  flag  FI  is  set  and  the  next  closest 
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value  of  VCV  Is  selected.  If  the  selected  value  Is  the  first  or  last  in 
the  array,  the  next  value  down  or  up  the  list  is  selected.  This  allows 
calculation  of  dC/dV. 

Line  13700-15500:  Search  of  the  array  CCV(I).  Beginning  with 
the  Jth  value  selected  intthe  previous  block,  values  of  CCV  are  identified 
which  are  at  least  1%  but  not  more  than  10%  different  from  the  Jth 
value.  If  the  search  does  not  find  values  to  fit  this  criteria,  flag  F 
is  set  and  first  CCV  value  is  selected  which  is  at  least  1%  different 
from  the  Jth  value.  These  are  designated  CCV(U)  and  CCV(L) . 

Line  15600-15900:  Values  of  dC/dV  (DCDV)  and  fixed  charge  con¬ 
centration  (NB(I-l) )  are  calculated.  Note  that  during  each  iteration, 
the  1-1  value  of  NB  and  NT  are  calculated. 

Line  16000-16300:  Debugging  print  statements  for  NB  and 
related  variables.  These  lines  may  be  deleted  if  desired. 

Line  16400:  In  the  case  1*1,  skip  the  blocks  concerning  trap 
concentration  calculation. 

Line  16500-19700:  Calculate  the  location  of  the  Fermi  level. 
This  is  accomplished  by  solving  the  transcendental  equation: 

N  exp(-(E  -E  )/kT)  *N./l+2exp(  (E.-E.)/kT)  +N  exp(  (E  -E.)/kT)  A.l 
c  cv  a  ra  v  vi 

The  value  of  E^  is  systematically  varied  until  the  difference  between  left 
and  right  hand  sides  of  equation  A.l  is  less  than  0.01%  of  E^,  of  until 
50  Iterations  are  completed.  In  the  latter  case,  an  error  message 
(line  19200)  is  printed  and  the  program  continues. 

Line  19800-21800:  This  block  calculates  the  depletion  width 
W(I) .  A  basic  derivation  of  the  expression  used  is  given  here.  We  begin 


with  the  egression  for  depletion  width  for  e  Schottky  barrier: 

-L-2-  (V  -  ~  -  V)  -  /"  x  N  (x)  dx  A.  2 

q  D1  H  O  D 

where 

Vbi  '  «bn  +i*  -  (Ec-Ef)/<l  A- 3 

Here  is  the  built-in  voltage,  V  is  the  applied  bias,  <{>bn  is  the 
Schottky  barrier  height,  A$  is  the  image  force  barrier  lowering,  N,  is 

D 

the  background  fixed  charge  concentration.  We  may  generalize  this 
equation  to: 

~TS~  <vbi  -  T  ‘  V  '  *  "i  W  A-4 

i-1 

£iit  &  -*■  « 

where  N^(Vi)  is  the  background  fixed  charge  concentration  at  the  deple¬ 
tion  depth  corresponding  to  V..  If  we  let  (e  e  )/q  *  b  and  V.  .  -  (kT/q) 

l  r  o  dj. 

Vj  “  vn’  we  wa^  write: 


bvx  -  Nb(Vx)  Wx 


bv2  "  Nb(Vl}  W1  +  Nb(V2)  [W2  ”  W2W1] 


bv  *  N,  (V.)W?  + . +  N,  (v  )  [W2  -  W  W  .] 

noil  bn  n  nn-i 

Substracting  the  last  two  lines  we  have: 

bVbVi- VV  ‘“n-VA1  *-6 

which  may  be  solved  quadratically  for  if  and  Nb(Vn>  are  known. 

W,  is  given  by  the  first  of  equations  A. 5.  This  scheme  is  implemented  in 


I 


--  w..  V.  .* 


103 


lines  20100  through  20600.  is  calculated  in  line  20000.  In  order 

to  obtain  the  best  estimate  of  W,  A$  is  calculated  iteratively  (DO  loop 

line  19900)  based  upon  the  latest  value  of  U  using: 

n 


A*±  - 


lE  “ 


4ire  e 
o  r 


,  2  2 
'■4ire  e 
o  r 


[E  < V V Wj  +  V V  <“| -Wj  _ J.)  )  ] 


i 


•  •  * 

V; 


p 


New  values  of  are  found  based  upon  this  value  of  A<|>^.  This  processes 

continues  until  the  difference  between  successive  calculations  of  W  is 

n 

.9 

less  than  10  cm.  We  assign  the  term  in  brackets  the  label  ACCDPT.  The 
value  of  ACCDPT  is  accumulated  for  the  next  profile  step  In  the  variable 
ADDCP.  Print  statements  21300-21600  were  for  debugging  and  are  not 
active. 

Line  21900-22300:  Calculate  average  value  of  & (AC/C)  (DELC(I-l). 

The  1-2,  1-1  and  Ith  values  of  DELC  and  VDELC  are  used.  If  I»2  or 
CPTS+1,  only  the  1-1  and  CPTS  values  of  DELC  and  VDELC  respectively  are 
used. 

Line  22400:  Calculate  trap  concentration.  See  section  2. 3. 2. 2. 

Line  22500-23500:  Print  results  of  this  profile  step. 

Line  23600-27700:  Termination  steps.  The  program  asks  if  add¬ 
itional  profiles  or  other  devices  are  to  be  analyzed.  If  so,  execution 
returns  to  lines  7000  and  1800  respectively. 

This  program  was  specifically  written  for  silicon  Schottky  barriers. 
Modification  for  other  semiconductor  types  requires  (1)  replacement  of 
line  9700  with  an  expression  to  calculate  E^  as  a  function  of  temperature 


00100 

C 

TRAP3I. 

00200 

c 

00300 

c 

00400 

c 

THIS  PROGRAM  COMPUTES  TRAP  DEPTH  PROFILES  PROM  A 

00500 

c 

SERIES  OP  DLTS  SCANS  OP  VARYING  BIAS  REDUCTING 

00600 

c 

PULSE  HEIGHTS.  PROGRAM  IS  DESIGNED  POR  3.B.  ON 

00700 

c 

N-TYPE  MATERIAL.  PROGRAM  COMPUTES  DEPLEATION 

00800 

c 

WIDTH  PROM  FIXED  CAPACITANCE  DATA  WHICH  YOU 

00900 

c 

MUST  PROVIDE.  THIS  VERSION  READS  C(V)  DATA  PROM  A  DISK 

01000 

REAL  S (50 ) , CVF( 50 ) , DISP ( 50 ) 

01100 

REAL  DELC ( 50 ) , DELCT (50 ) , VB( 50 ) , W(0 : 50 ) 

01200 

REAL  AREA , ER , MDSE , MDEH , EGEP , VBIAS , CBIAS , VR 

01300 

REAL  MI,N3,NV,EG,EF, EPDIPP , INC , EDONN , DELPHI , DPHIBN, WOLD 

01 400 

REAL  ACCDP.DCDV, CTEMP.PTI ,  FT  2,  PT 3, ACCDPT 

01500 

REAL  A(3O),NT(50),NB(0:50),CCV(lO0),VCV(lOO),VDELq(5O) 

01600 

INTEGER  CFTS,F,F1 , J, K, I,U, L, Q,X,N, F2, NP0INT,0,T 

01700 

DOUBLE  PRECISION  B 

01300 

1 

WRITE(5,10) 

01900 

10 

FORMAT (IX, 'ENTER  ONE  LINE  COMMENT  ',/) 

02000 

READ(5, 20)  A 

02100 

20 

FORMAT (80A1) 

02200 

PRINT  30, A 

02300 

30 

FORMAT (1X,80A1) 

02400 

WRITE (5, 40) 

02500 

40 

F0RMAT(1X, 'ENTER  DEVICE  AREA  (CM**2)  -  ',*) 

02600 

READ(5, 50)  AREA 

02700 

50 

FORMAT (G) 

02800 

WRITE(5,60) 

02900 

60 

FORMAT (IX, 'ENTER  RELATIVE  DIELECTRIC  CONST.  -  ',$) 

03000 

READ(5,50)  ER 

03100 

VRITE(5,65) 

03200 

65 

FORMAT (IX, 'ENTER  DELTA  PHI  B-N  (VOLTS )-',*) 

03300 

READ(5,50)  DPHIBN 

03400 

WRITE (5, 70) 

03500 

70 

FORMAT (IX, 'DENSITY  OF  STATES  EFFECTIVE  MASSES  OF  COMMON 

03600 

C 

'  SEMICONDUCTORS:'/ 

03700 

C 

'  ELECTRONS  HOLES'/ 

03800 

c 

'SI  1.1  .557 

03900 

c 

'  GAAS  .068  .  57 

04000 

c 

'  GAASP  40 %  .039  .5'/ 

04100 

c 

'  GAASP  INDIRECT  1.20  .57 

04200 

c 

'GAP  1.20  .57) 

04300 

WRITE (5. 80) 

04400 

80 

FORMAT (1X,'MDS  FOR  ELECTRONS*' ,$) 

04500 

READ(5,50)  MDEE 

04600 

'WRITE  (5, 90) 

04700 

90 

FORMAT (IX.'MDB  FOR  HOLES-', *) 

04300 

READ(5, 50)  MDEH 

04900 

WRITE (5, 91  ) 

05000 

91 

F0RMAT(1X, 'DONOR  ACTIVATION  ENERGY  (EV )-',*) 

05100 

READ(5,50)  EDON 

05200 

WRITE(5, 100) 

05300 

100 

FORMAT (IX, 'ENTER  NAME  OF  C(V)  DATA  FILE:  ',$) 

05400 

READ(5,105)  B 

FORMAT (A9) 

OPEN  (UNIT-20, DEVICE-'DSK' , ACCESS-'SEQIN' , FILE-B) 

DO  130  1-1,100 
READ(20,110)  CTEMP.VCV(I) 

FORMAT (20) 

IF  (CTfflP.EQ.O. AND.VCV(l) .EQ.O)  120,125 

NPOINT-I-1 

GO  TO  138 

CCV(l)-CTEMP-0.68 

CONTINUE 

WRITE(5,131) 

FORMAT(lX,'C(V)  DATA  FILE  TOO  URGE  —  CONTINUING  ANYWAY',/) 
GO  TO  140 
WRITS(5, 139) 

FORMAT  (1X,'C(V)  DATA  READ  A-OK',/) 

WRITE(5, 141 ) 

FORMAT(l X, 'ENTER  1  LINE  COMMENT',/) 

RSAD(5,20)  A 
PRINT  30,  A 
PRINT  145 

FORMAT (2X,'F',3X,'V-REV',8X,'S.S.  CAP' ,10X, 'EC-SF' , 

C  12X, 'DELPHI', 7X,'AVG.  D(DC/C)/DV',4X,'DEP.  WIDTH', 

C  5X, 'DOPANT  CONCNT.', 5X, 'TRAP  C0NCNT.',/,5X, '(VOLTS)', 

C  9X, '(PF)',12X,'(EV)',13X, '(VOLTS)', 9X,'(1 /VOLTS)', 

C  9X,'(CM)',12X,'( CM**-3 ) ' , 1 0X , ' ( CM**-3 ) ' ) 

WRITE(5,150) 

F0RMAT(1X, 'VALUE  OF  DC  BIAS  (VOLTS)-', $) 

READ(5,50)  VBIAS 

DO  160  1-1  .NPOINT 

IF  (VBIAS.EQ.VCV(I))  GO  TO  180 

CONTINUE 

WRITE(5, 170) 

F0RMAT(1X,'VBIAS-VCV(I)  NOT  FOUND.  C(VBIAS)  SET  TO  C(2)',/) 

PRINT  170 

CBIAS«CCV(2) 

GO  TO  181 
CBIAS-CCV(I) 

WRITS(5,190) 

FORMAT( IX, 'ENTER  TEMP  OF  PEAK(DEC  K)  ',$) 

READ(5,50)  TIMP 

NC-4.8290E15*MDEE**1 .5*TEMP**1.5 
NV-4.8290E15*MEGH**1 ,5*TEHP**1 .5 
EG-1  .16-(7.02E-4*TEMP**2)/(1 108+TEMP) 

WRITE(5,210) 

F0RMAT(1X, 'ENTER  DELTA  CAP  DATA.  4  OF  PTS-',S) 

READ(5,50)  CPTS 

W0LDM3.0O1E-5 

ACCDPKD 

ACCDPT-0.0 

W(0>0 

DO  510  1-1  ,cprs*i 

F-0 

F1-0 

IF  (I«EQ.CFTS*1  )  GO  TO  245 
WRITE(5, 220)  I 
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M 

11000 

11100 

220 

B 

11200 

11300 

11400 

230 

11500 

11600 

1 1700 
11800 

240 

I*  , 

11900 

235 

12000 

12100 

12200 

*% 

• 

& 

12300 

12400 

12500 

12600 

245 

I* .  • 

12700 

12300 

250 

12900 

e 

13000 

13100 

13200 

13300 

253 

1  * 

13400 

■ 

13500 

13600 

13700 

255 

p 

13800 

260 

13900 

14000 

?:;• 

14100 

270 

c.> 

..  «l 

m 

14200 

14300 

14400 

280 

K 

14500 

290 

14600 

300 

14700 

310 

h  7 

u . 

14800 

14900 

Lm* 

15000 

320 

15100 

330 

15200 

15300 

15400 

340 

15500 

350 

15600 

360 

15700 

15800 

iZj 

15900 

'•I'-' 

16000 

362 

16100 

16200 

361 

d 

16300 

16400 

F0RMAT(1X,'BIAS  REDUCTING  PULSE(V0LTS)  V(', 13, ')-',$) 
READ(5,50)  VR 
VDELC(l)-VBIAS+VR 
WRITS(5,230)  I 

F0RMAT(1X, 'ENTER  VAL  #',I3,'  OF  DISP  (CM),  SENS.  (V), 

C  CONVSR.  FACT.  (V/PF)',/) 

READ(5,240)  DISP(l),S(l),CVF(l) 

FORMAT (3(0)) 

WRITE(5,235)  I , VR, DISP(l) ,S(l) ,CVF(I) 

F0RMAT(lX, 'VALUE  #  ',13, 'OF  PULSE-', FI 2. 5,/,'  OF  DISP-', 

C  F12.5 ,/,'  OF  SENS-', FI 2. 5,/,'  OF  CON.  FACT.-',F12.5,/, 

C  '  OK?  —  1 -YES,0-N0:  ',*) 

READ  (5,50)  T 
IF  (T.EQ.O)  GO  TO  219 

DELCT(I)-(S(I)*DISP(I)*.5650*.5)/(CVF(I)*CBIAS) 

IF  (I.EQ.1)  VDELC ( 0 )-VBIAS 
DO  250  J-1.NP0INT 

IF  (ABS (VDELC (1-1  )-VCV(J)) .LE. IE-5)  GO  TO  260 

CONTINUE 

FI-1 

MI -VDELC  ( 1-1  ) 

DO  255  J-1 , NPOINT 

IF  (ABS (VDELC (1-1  )-VCV(j)) .LT.MI)  253,255 
MI -ABS  (VDELC  (1-1  )-VCV(j)) 

Q-J 

CONTINUE 

IF  (Q.EQ. NPOINT)  J-Q-1 
IF  (Q.EQ.  1 )  J-2 
DO  270  K-1, NPOINT 
IF  (K+J.EQ. NPOINT)  GO  TO  280 

IF  (AB3(CCV(K+J)-CCV(J)).<7T.(0.01*CCV(J)))  GO  TO  290 

CONTINUE 

U-KfJ 

F-1 

GO  TO  310 

IF  (ABS(CCV(K*J)-CCV(j)).LE.(0.1*CCV(j)))  300,230 
U-K+J 

DO  320  K-1 .NPOINT 
IF  (J-K.EQ.1 )  GO  TO  330 

IP  (ABS(CCV(J-K)-CCV(J)).CT.(0.01#CCV(J)))  GO  TO  340 

CONTINUE 

L-J-K 

F-1 

GO  TO  360 

IF  (ABS(CCV(J-K)-CCV(J)).LE.(0.1*CCV(J)))  350,330 
L-J-K 

dcdv-((abs(ccv(u)-ccv(j)))/(  ABS  (vcv(u)-vcv(j  )))•*• 

C  (ABS(CCV(j)-CCV(L)))/(ABS(VCV(j)-VCV(L))))/2 
NB(I-1  )-CCV(j  )##3/(  1.41 8E-8*ER*AREA**2*DCDV) 

IF  (NB(l-1  ).LT.O)  NB(I-1  )-ABS(NB(l-1  )) 

PRINT  361,NB(I-1  ) ,DCDV,CCV(J)  ,CCV(U)  ,CCV(D 
WRITE  (5,361)  NB(I-1 ),DCDV,CCV(J),CCV(U),CCV(L) 

FORMAT('  NB(I-1  )-' ,1PE1 1 .3, '  DCDV-',1PE1 1 .3,'  CCV(J)-', 

C  2PE11.3,'  CCV(U)-',2PE11.3,'  CCV(L)-' ,2PE1 1 .3) 

IF  (I.EQ.1)  GO  TO  510 
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m 


16500 

16600 

16700 

16800 

16900 

17000 

17100 

17200 

17300 

17400 

17500 

17600 

17700 

17800 

17900 

18000 

18100 

18200 

18300 

18400 

18500 

18600 

18700 

18800 

18900 

19000 

19100 

19200 

19300 

19400 

19500 

19600 

19700 

19800 

19900 

20000 

20100 

20200 

20300 

20400 

20500 

20600 

20700 

20300 

20900 

21000 

21100 

21200 

21300 

21400 

21500 

21600 

21700 

21300 

21900 


INC-0.10«EG 

F2-1 

EFO.50*EG 
DO  390  0-1 ,50 

PEI -EXP  (-1  .0*(EG-EF)/(8.6l71E-5«TBtP)) 

PT2-NB(l-1  )/NC 

PT3-(1  -0/(1 .  0*2. 0*(EXP((EF-EG«*D0N)/(8. 6171  B-5*T»P)  )  )  )  ) 

EPDIPP-PT  1-PT2*PT3 

IP  (0.EQ.1)  GO  TO  355 

GO  TO  365 

SIGN-EFDIFF 

GO  TO  385 

IP  (P2.EQ.1)  370,380 

IP  ( (SIGN. LT. 0. AND. EPLIPP . GT . 0) . OR. (SIGN. OT . 0. AND. EPDIf F. 

LT. 0 ) )  371,372 
F2--1 

SIGN-EFDIFP 
INC-O. 1*INC 
GO  TO  372 

IP  ( (SIGN. LT. 0. AND. EFDIFF. GT . 0) .OR. ( SIGN. GT . 0. AND. 

EFDIFF.LT. 0))  381,372 
F2-1 

SICN-EFDIFF 
INCK).  1  *INC 

IF  ( ABS (EFDIFF) .LT. (EG*1  E-4 ) )  400,385 
EF-EF+(INC*F2) 

CONTINUE 

WRITE(5,395)  I, EFDIFF, EG, EF,NC,Pri,Pr2,Pr3,NB(l-1  ),DCDV 
FORMAT(lX,'EF  ROUTINE  UNABLE  TO  CONVERGE  IN  50  ITERATIONS. 
',/,1X,'QUIT  AT  POINT:',  13, 2X, 'EFDIFF-' ,E1 1 -3,2X, 'EG-' ,E1 1 .3, 
'  EF-' ,E1 1.3,'  NC-' ,E1 1.3,'  PT1-',E11.3,'  PT2-',E11.3, 

'  PT3-',E1 1 .3, '  NB(l-1  )-' ,E1 1.3,'  DCD-'  ,E1 1.3) 

PRINT  395,I,EFDIFF,EG,EF,NC,Pri,Pr2,Pr3,NB(l-1  ),IJCDV 
IF  (I.EQ.2)  DELPHI-0.01 
DO  450  0-1,50 

VB(  I  -1  )-  DELPH I +DPH IBN  -  ( EG-EF)  -VDELC  ( I  -1  )-TENP*8.6l  71 B-5 
IF  (I.EQ.2)  410,420 

W(I-1  )-SQRT(ER*5.527E5*V3(l-1  )/NB(l-1  )) 

GO  TO  440 

TERN-(SQRT(W(l-2)**2+4.0#5.527E5*ER*(VB(l-1 )-VB(l-2)) 

/NB(I-1  ))/2.0) 

V(l-1  )-W(l-2)/2+TERM 
IF  (W(l-1  ).GE.O)  440,430 
W(I-1  )*W(l-2)/2-TERM 

IF  (ABS(W(I-1  )-WOLD).LT.0.0001E-5)  CO  TO  460 
WOLD-W  ( I  -1  ) 

ACCDPT "(W(  1-1  )-W(l-2  ))*NB(l-1  )+ACCDP 
DELPHI-4 . 325E-8*SQRT( ACCDPT ) 

WRITE  (5,441)  W(  1-1  ),VB(I-1  ), TERN,  ACCDPT,  DKLPI  l ,  VDELC(  I -1  ) 
PRINT  441 ,  W(l-1  ) , VB(l-1  ),TERM,  ACCDPT, DELPHI, VDBLC(l-1  ) 

FORMAT  ('  W(l-1  )-',E12.4,'  VB(l-1  )-'  ,E1 2. 4, '  TERM-' ,E12. 4, 

'  ACCDPT-', El  2. 4,'  DELPHI-', E12. 4,'  VDELC (1-1 )-' ,E12.4) 

CONTINUE 

ACCDP-ACCDPT 

I?  (I.EQ.2.0R. I.EQ.CPTS+1  )  470,430 


22000 

470 

DBLCd-1  )»(DELCT(l)-DELCT(l-1  ))/(VDELC(l)-VDELC(l-1  )) 

22100 

GO  TO  490 

22200 

480 

DELC(I-1  )-(((delct(i)-delct(i-i  ))/(VDELC(l)-VDELC(l-1  ) 

22300 

C 

))*((DELCT(I-1  )-DELCT(l-2))/(VDELC(l-1  )-VDELC(I-2))))/2 

22400 

490 

ST (1-1  )-ABS(l  .809E-6«NB(l-1  )*NB(0)*W(I-1  )**2*DELC(I-1  )/ER) 

22500 

IF  (F.NE. 1 .AND. FI .NE. 1 )  X-0 

22600 

IF  (F.NE.1.AND.F1.EQ.1)  X-1 

22700 

IF  (F.EQ.1.AND.F1.NE.1)  X-2 

22800 

IF  (F.EQ.1.AND.F1.EQ.1)  X-3 

22900 

BGEF-EG-EF 

23000 

WRITE  (5,500)  X,VDELC(l-1  ) ,CCV(j) ,EGE7, DELPHI, DELC(I-1  ), 

23100 

c 

V(I-1  ),NB(l-1  ),NT(I-1  ) 

23200 

500 

F0RMAT(2X,I1,2X,F7.3,7(5X,1PE12.5)) 

23300 

FRIST  500,X,VDELC(I-1  ),CCV(J) ,EGEF, DELPHI, DEM (1-1  ) , 

23400 

c 

V(I-1  ),SB(I-1  ),NT(I-1  ) 

23500 

510 

CONTINUE 

23600 

PRINT  511, AREA 

23700 

5 11 

FORMAT (///, IX, 'DEVICE  AREA-" , 1PE1 2. 5, 'CM**2 ') 

23800 

PRINT  512.ER 

23900 

512 

FORMAT (l  X, 'RELATIVE  DIELECTRIC  CONST. «' ,F6. 2) 

24000 

PRINT  513, DPHIBN 

24100 

513 

FORMAT (l X, 'DELA  PHI  B-N-' ,F6. 2, '(VOLTS)') 

24200 

PRINT  514,MDEE,MDEH 

24300 

514 

F0RMAT(1X, 'DENSITY  OF  STATES  EFFECTIVE  MASS  FOR  ELECTRONS-' 

24400 

c 

F4.2,'  AND  OF  HOLES-' ,F4. 2) 

24500 

PRINT  515.ED0N 

24600 

515 

F0RMAT(1X, 'DONOR  ACTIVATION  ENERGY-', F5. 3, '(EV)') 

24700 

HUNT  51 6,  TEMP 

24800 

516 

F0RMAT(1X, 'TEMPERATURE  OF  PBAK-',F6.2,'(DEG  K)') 

24900 

PRINT  517.VBIAS 

25000 

517 

FORMAT (l X,'DC  REVERSE  BIAS-' ,F7. 2, '(VOLTS)') 

25100 

PRINT  520 

25200 

520 

FORMAT(/,'  FUG  CODES:  0  -  NO  DIFFICULTY;  1  -  COULD 

25300 

c 

NOT  MATCH  A  C(V)  DATA  POINT  TO  THIS  VALUE  OF  VREV; 

25400 

c 

2  -  DELTA  C/DELTA  V  CRITEIA  NOT  MET;',/,'  3  -  BOTH 

25500 

c 

CONDITIONS  1  A  2  APPLY.' ,///,30X, 'CAPACITANCE  DATA 

25600 

c 

US  ED  ABOVE ' ,/ , 6X , 'V-REV' , 8X , 'HEIGHT' , 9X , 'SENSITIVI TY' , 

25700 

c 

5X, 'CONVER.  FACT. ' ,8X , 'DELC/C' ,/ , 5X, ' (VOLTS' ,9X , 

25800 

c 

'(CM) ' ,  1 1 X , '  (VO  LTS ) ' ,  1 1 X , '  (V  /PF ) ' ,  1 1 X , '  (PF/V ) ' ) 

25900 

DO  540  M-1  ,CPTS 

26000 

PRINT  530,VDBLC(M),DISP(M),S(M),CVF(M),DBLCT(M)  . 

26100 

530 

F0RMAT(5X , F7. 3, 5 (5X , 1 PE  1 2 . 5 ) ) 

26200 

540 

CONTINUE 

26300 

PRINT  545 

26400 

545 

FORMAT(//////) 

26500 

WRITE(5,550) 

26600 

550 

F0RMAT(l X,'IS  THERE  ANOTHER  TRAP  PEAK  FOR  THIS  DEVICE?',/, 

26700 

c 

1 X, '1 -YES,0-N0' ,/) 

26800 

READ(5,50)  S 

26900 

IF  (N.EQ.1)  GO  TO  140 

27000 

VRITE(5,560) 

27100 

560 

70RMAT(1X,'D0  YOU  WISH  TO  ANALYSE  ANOTHER  DEVICE?',/, 

27200 

c 

IX, '1 -YES,  O-NO',/) 

27300 

READ (5, 50)  N 

27400 

I?  (N.EQ.1)  570,580 
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